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ABSTRACT 

i 

i 

V  . 

A  new  very  single  and  dependable  field  sediment  sampler  has 
been  developed.  It  is  a  multi -siphon  device  that  samples  suspended 
sand  from  seven  elevations  simultaneously.  The  65  concentration 
profiles  that  have  been  obtained  so  far  are  presented,  complete 
with  hydrodynamics  and  bed  sediment  data.  The  profiles  were 
measured  under  breaking  and  non -breaking  waves  in  a  wide  range 
of  coastal  environments  all  over  Australia.  These  data  together 
with  a  review  of  oscillatory  boundary  layer  flow  and  the  motion  of 
suspended  sand  in  an  accelerated  non  uniform  flow  lead  to  a  range 
of  new  insights  into  the  processes  of  sediment  entrainment  under 
vaves . 
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1.  INTRODUCTION 

Sediment  transport  in  the  coastal  environment  is  often  a 
result  of  the  combined  actions  of  waves  and  currents.  A  steady 
current  must  at  first  sight  seem  more  efficient  as  a  sand  transporter 
than  the  back  and  forth  motion  of  waves;  but  the  process  of  picking 
the  sand  up  from  the  bed  is  dominated  by  the  waves.  The  efficiency 
of  waves  with  respect  of  entrainment  of  sediment  is  due  to  two  main 
characteristics  of  oscillatory  water  motion. 

Firstly,  the  boundary  layer  of  an  oscillatory  flow  is  thin: 
a  few  centimetres  for  typical  wave  periods,  which  means  that  wave 
induced  near  bed  velocities  and  bed  shear  stresses  become  large. 

Secondly,  the  periodic  flow  reversals  cause  the  boundary 
layer  to  separate  from  the  bed  twice  every  period  and  send  dense 
clouds  of  suspended  sediment  upwards.  This  is  the  first  step  of  the 
entrainment  process  by  which  the  suspended  sand  reaches  elevations 
comparable  to  the  wave  boundary  layer  thickness.  From  this  level 
it  may  either  start  to  settle  out  again  or  rise  further  due  to 
the  mixing  action  of  currents  or  wave  breaking. 

The  present'  report  gives  a  detailed  discussion  of  these  pro¬ 
cesses  based  on  the  literature  and  previous  laboratory  studies  as 
well  as  a  large  amount  of  field  data  collected  by  the  Coastal 
Studies  Unit  over  the  last  two  and  a  half  years. 

Chapter  2  deals  with  the  hydrodynamic  problems  involved  in 
sediment  transport  calculations,  mainly  the  structure  of  wave 
boundary  layers.  Special  emphasis  is  put  onto  developing  formulae 
relating  quantities  such  as  boundary  layer  thickness,  bed  shear 


stress  and  the  turbulence  structure  with  very  few  assumptions  about 
the  nature  of  the  flow. 

Previous  experimental  studies  fall  into  three  categories  and 
the  linking  together  of  their  results  has  so  far  been  in conflate. 
Detailed  velocity  measurements  like  those  of  Kalkanis  (1964)  and 
Jonsson  and  Carlsen  (1976)  give  good  information  about  the  average 
character  of  the  flow  and  the  shear  stresses,  but  they  get  no  direct 
information  about  the  turbulence  structure.  Such  information  has 
become  available  later  through  the  measurements  of  MacDonald 
(1977) ,  Nakato  et.  al.  (1977) ,  Du  Toit  and  Sleath  (1981)  and 
Kemp  and  Simons  (1982) .  However,  all  of  these  except  from  Nakato 
et.  al.  used  artificial  solid  beds.  Natural  sand  beds  were  used 
by  Cars tens  et.  al.  (1969)  and  Dofquist  (1980)  who  measured  energy 
dissipation  and  bed  friction.  The  major  aim  of  the  analysis  is  to 
make  the  latter  friction  measurements  useful  for  prediction  of  eddy 
viscosities  and  vertical  length  scales  over  natural  sand  beds. 

Chapter  3  is  an  analysis  of  the  motion  of  suspended  sediment 
particles  in  an  accelerated,  non  uniform  flow  field.  The  approach 

to  the  solution  is  new  in  that  it  splits  the  sediment  particle 

->•  ^  ^  ^  ^ 
velocity  u  into  three  parts  from  the  outset,  u  ■  u+v+w.  u  xs  the 
P  P 

water  velocity,  w  is  the  still  water  settling  velocity  and  v  is  the 
extra  velocity  which  results  from  the  fluid  accelerations.  Special 
attention  is  paid  to  the  case  of  sand  in  vortices  because  this  is 
the  key  to  understanding  the  entrainment  mechanism  over  rippled  beds 
under  waves,  and  the  resulting  sediment  distributions.  The  simplest 
approximation  u^  *  u+w  and  the  very  simple  kinematic  arguments  are 
used  to  show  the  essence  of  the  trapping  mechanism  which  makes  it 
possible  for  grains  of  nearly  all  sizes  to  travel  together  inside 
migrating  vortices.  A  process  which  leads  to  very  similar  dis- 
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tributions  of  all  grain  sizes  over  ripples.  This  process  is  funda¬ 
mentally  unlike  gradient  diffusion  which  leads  to  vastly  different 
distributions  for  different  grain  sizes.  The  fact  that  sand  and 
vortices  travel  together  also  explains  why  the  distributions  of  sand 
and  turbulence  are  very  similar  as  observed  by  Nakato  et.  al.  (1977) 
and  MacDonald  (1977) . 

Chapter  4  describes  the  experimental  equipment  and  experimental 
technique.  The  sediment  sampler  was  designed  to  sample  a  complete 
profile  of  time  averaged  concentrations  simultaneously.  This  was 
done  in  order  to  overcome  the  large  scatter  problem  encountered  in 
previous  field  studies,  e.g.  Fairchild  (1977)  and  Kana  (1979). 

After  obtaining  well  defined  concentration  profiles  the  problem  is 
to  extract  parameters  that  are  both  practically  and  theoretically 
well  defined,  which  can  in  turn  be  related  to  the  flow  of  sediment 
characteristics.  In  the  present  study  we  extract  these  parameters 
from  the  profiles:  C0,  eQ  and  Eg.  Cq  is  the  concentration  at  the 
bed  level.  Eg  represents  the  diffusivity  inside  the  wave  boundary 
layer  and  sg  represents  the  larger  diffusivity  outside  the  boundary 
layer,  due  to  wave  breaking  or  currents.  Some  concentration  profiles 
show  an  apparent  increase  in  diffusivity  or  vertical  length  scale 
with  distance  from  the  bed,  although  there  is  neither  breaking  waves 
nor  steady  currents  around  to  produce  significant  turbulence  outside 
the  boundary  layer.  This  apparent  increase  is  due  to  the  variability 
of  settling  velocity  in  the  suspended  material.  A  method  for 
correcting  the  observed  diffusivities  for  this  effect  is  devised. 

It  is  based  on  the  diffusion  equation  and  uses  a  continuous  settling 
velocity  distribution  which  leads  to  vastly  simpler  calculations 
than  previous  models  that  work  with  discrete  size  classes. 


Chapter  5  is  a  heuristic  discussion  and  analysis  of  the  field 
measurements  of  the  present  study  together  with  the  results  of 
some  previous  laboratory  investigations.  The  clearest;  results 
are  obtained  in  relation  to  the  near  bed  sediment  distributions. 

Over  fully  developed  vortex  ripples  we  find  exponential  con¬ 
centration  profiles 

c(z)  ■  C  exp  (-z/£)  Cl.l) 

o  s 

where  the  length  scale  Ag  is  approximately  equal  to  the  ripple  height 
Different  grain  sizes  follow  very  similar  distributions  because  the 
entrainment  mechanism  is  convective  as  described  in  Chapter  3.  The 
sediment  gats  trapped  in  the  lee  vortices  while  these  are  being 
formed,  and  when  the  vortices  are  released  and  travel  upwards,  they 
act  as  elevators  for  all  grain  sizes  alike. 

Over  flat  beds  and  round  crested  bed  forms  the  vertical  length 
scale  of  the  concentration  profiles  is  typically  ten  times  smaller 
than  over  ripples,  so  it  is  rather  difficult  to  obtain  detailed 
information  about  the  profile  structure.  However, the  present  data 
give  a  good  indication  that  the  lower  part  of  the  profiles  scale  on 
the  hydraulic  roughness  length  of  the  flat  bed,  and  that  the 
entrainment  mechanism  within  the  wave  boundary  layer  is  predominantly 
convective  because  the  observed  near  bed  length  scales 
£g  ("Cg/w)  depend  little  on  the  settling  velocity. 

The  concentration  magnitudes,  CQ,  is  somewhat  more  difficult 
to  predict.  It  is  probably  valid  to  assume  that  CQ  depends  mainly 
on  the  skin  friction  Shields  parameter,  9*.  However,  a  proper  cal¬ 
culation  of  9'  for  the  laboratory  data  is  difficult  because  most  of 
them  are  in  the  transition  region  between  laminar  and  turbulent  flow. 


6 


2. 


HYDRODYNAMICS 


The  water  motion  in  the  nearshore  area  and  the  3urfzone  involves 
velocities  that  vary  over  a  large  range  o£  time  scales,  from  diurnal 
tides  to  harmonics  of  wind  waves  with  periods  less  than  a  second. 

It  is  often  convenient  to  split  the  velocity,  u,  in  two  parts 

u  »  u  +  u  (2.1) 

where  u  is  the  steady  or  quasi  steady  part  and  u  is  periodic  with 
zero  mean. 

The  cutoff  period,  above  which  everything  is  considered  quasi 
steady  depends  on  the  problem  under  consideration.  Sometimes  every¬ 
thing  above  wind  wave  periods  (a  10  seconds)  can  be  considered  steady 
in  other  cases  even  the  tide  must  be  considered  unsteady. 


WAVES 

The  most  important  part  of  the  water  motion  is  due  to  wind  waves 
and  swell  with  periods  between  5  and  20  seconds.  At  present  there  are 
no  theories  for  non  periodic  waves  so  for  any  theoretical  considerations, 
waves  are  assumed  periodic.  This  is  of  course  an  idealized  description 
since  the  wave  period  is  normally  quite  variable  just  like  the  wave 
height  H  in  the  surf  zone.  The  pressure  record  in  Figure  2.1  gives  an 
idea  of  this  variability.  It  is  necessary  however  for  any  theoretical 
work  to  represent  the  incoming  waves  by  a  monocromatic  wave  train  with 
constant  height.  With  respect  to  sediment  transport  the  most  relevant 


wave  period  is  the  one  determined  from  zero  crossings  of  the  near  bed 


Figure  2.1;  Instantaneous  surface  elevations  measured  by  a  pressure 
transducer.  The  waves  we  spilling  breakers  under  a  slight  onshore 
wind  (TEST  21,  Broome  1981) .  Note  that  the  time  progresses  from  right 


to  left.  With  typical  records  like  this,  the  meanings  of  such  terms 
as  wave  period  and  wave  height  are  no  longer  trivial. 
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velocity  because  these  flow  reversals  cause  a  break  down  of  the  boundary 
layer  structure. 

Individual  waves  are  best  characterized  by  their  zero  downcrossing 

height  H,  and  for  an  ensemble  of  waves  one  may  use  either  the  average 

H,  the  root-mean-square  H  or  the  significant  wave  height  H  ,  which 

ms  s 

is  the  average  of  the  highest  1/3  of  the  waves.  The  relevance  of  each 

of  these  measures  depends  on  the  phenomenon  described.  There  is  some 

evidence  that  H  is  superior  for  describing  ripple  geometry  (Nielsen 
s 

1981)  while  H  seems  to  be  the  best  for  describing  sediment  con— 
nns 

centrations  (Rasmussen  and  Fredsoe  1981) . 

Since  nearshore  waves  are  often  highly  non  linear  and  irregular 
the  only  way  to  get  reliable  information  on  the  water  motion  at  a 
given  spot  is  by  measuring  it  locally. 

The  best  way  of  measuring  the  velocity  field  is  by  using  an  array 
of  good  flow  meters .  However,  one  may  also  derive  the  characteristic 
wave  induced  velocities  from  pressure  measurements,  and  as  long  as 
our  main  concern  is  the  velocity  right  at  the  bed,  we  may  use  linear 
wave  theory.  The  relations  between  wave  height  and  the  pressure  and 
velocity  fluctuations  for  a  sinus  wave  are 


?(z) 


H  cosh  kz 
°^2  sinh  kh 


(2.2) 


0(z) 


TTH  cosh  kz 
T  sinh  kh 


(2.3) 


and  for  the  near  bed  velocity  amplitude  u(o)  which  will  be  referred  to 
in  the  following  as  aw  we  get 


T  sinh  kh 


(2.4) 


27T/L 


the  important  parameter  kh  which  is  the  local  wave  number,  k  * 
times  the  water  depth^ is  determined  by  the  dispersion  relation 


k  h  ■  kh  tanh  kh  (2.5) 

o 

This  implicit  relation  has  traditionally  been  solved  numerically 
and  the  results  tabulated,  however  it  is  accurate  enough  and  more 
practical  -to  use  the  following  explicit  formula 

kh  -  >/5Th  [l  +  frch  +  fa  OcQh)*]  (2.6) 

where  k  is  the  deep  water  wave  number 
o 

k  -  2ir/Ti  -  4ir^gT»)  (2.7) 

o  O' 

Equation  (2.6)  is  exact  in  the  limit  k^h  -*■  o  and  the  relative  error 

is  less  than  0.0045  for  k  h  <  2.5.  In  deeper  water  we  recommend 

o 

the  asymptotic  expansion . 

kh  »  kjjh  (1  +  2e“23^ h)  (2.8) 

The  two  formulae  (2.6)  and  (2.3)  cover  the  whole  range  0  <  k  h  <  «® 
with  an  error  of  less  than  0.45  percent.  Relative  errors  for  the 
two  formulae  and  for  a  few  simpler  formulae  as  well  are  shown  in 
Figure  2.2.  Details  about  the  derivation  of  equation  (2.5)  from 
the  linear  wave  dispersion  relation  are  given  by  Nielsen  (1982). 
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WAVE  BOUNDARY  LAYERS 

The  wave  boundary  layer  is  intuitively  defined  as  the  layer  close 
to  the  bottom,  where  the  wave  induced  water  motion  is  noticeably 
affected  by  the  boundary.  This  layer  is  normally  very  thin,  i.e. 
a  few  millimetres  over  a  smooth,  solid  bed  and  a  few  centimetres  over  a 
flat  bed  of  loose  sand.  Bed  forms  like  ripples  will  change  the 
structure  of  the  boundary  layer  by  introducing  strong  rythmic 
vortices  and  the  boundary  layer  over  sharp  crested  ripples  will 
extend  to  a  height  of  four  or  five  ripple  heights  or  up  to  about 
50  centimetres  under  field  conditions. 

It  is  a  typical  feature  of  oscillatory  boundary  layers  that  the 
velocity  close  to  the  boundary  and  the  resulting  shear  stress  t(t) 
are  ahead  of  the  free  flow,  u.  (t) .  For  laminar  flow  the  shear 
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stress  is  45  degrees  ahead  and  for  rough  turbulent  flow  it  is  between 
20  and  30  degrees  ahead. 

The  shear  stress  in  a  simple  harmonic,  laminar  flow  is  simple 
harmonic,  but  in  turbulent  flow  the  variation  with  time  is  much  more 
complicated.  When  the  ratio  between  the  bed  roughness,  r  and  the 
semi  excursion,  a  are  small,  the  variation  of  T  (t)  is  still  quite 
smooth  and  rather  like  a  simple  harmonic.  This  is  the  case  for  the 
measurements  of  Jonsson  and  Carlsen  (1976)  where  r/a  was  only 
0.00a,  see  Figure  2.3.  For  flow  over  fully  developed  ripples, 
the  ratio  r/a  is  of  the  order  of  magnitude  one,  and  the  flow  near 
the  bed  is  dominated  by  the  rythmic  formation  and  release  of  strong 
vortices.  Lofquist  (1980)  measured  instantaneous  values  of  T  (t) 
under  such  conditions.  Figure  2.4  shows  s ome  of  his  results  and  we 
see  that  the  behaviour  is  completely  different  from  that  of  a  simple 
harmonic  or  from  that  of  sin  (cut  -  P)  |  sin  (cut  -  i)  j  which  has 
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Figure  2.3:  Tina  variations  of  the  bed  shear  stress  T(0,t) 
for  rough  turbulent  flow  over  relatively  small  roughness  elements 
r/a  a  .008.  From  Jonsson  and  Carlsen. 
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been  assumed  in  several  "theoretical"  studies. 

Another  typical  feature  of  oscillatory  boundary  layers  is  the 
"overshoot”  near  the  bed,  which  comes  about  because  the  velocity 
deficit  3  (t)  -  u(z,t)  has  the  nature  of  a  damped  wave  which 

altematingly  adds  to  and  subtracts  from  the  total  amplitude  u(z)  . 

See  Figure  2.5. 

The  above  mentioned  main  periodic  properties  of  the  oscillatory 
boundary  layers  have  been  studied  in  great  detail  by  Kalkanis  (1957 
and  1964) ,  Kamphuis  (1975)  and  Jonsson  and  Carlsen  (1976)  on  the 
experimental  side  and  by  Kajkara  (1968)  and  Jonsson  (1980)  on  the 
theoretical  or  speculative  side.  Their  work  is  useful  for  sediment 
transport  studies  because  they  provide  predictive  formulae  for  bed 
shear  stresses  and  boundary  layer  thicknesses.  Those  early  experimental 
studies  were  not  however  concerned  with  the  details  of  the  turbulence 
or  the  vortex  motion.  Nevertheless,  the  latter  are  very  important 
for  the  understanding  of  sediment  entrainment  processes.  The 
formation  and  strength  of  lee  vortices  over  large  ripples  was  studied 
and  described  by  Tunsdall  and  Inman  (1975) ,  who  found  that  the 
vortices  formed  twice  every  period  just  about  the  time  of  velocity 
extremes  in  either  direction  and  that  the  order  of  magnitude  of 
the  maximum  velocities  in  the  vortices  was  aw. 

The  turbulence  (or  velocity  fluctuation)  intensities  over  ripples 
in  purely  oscillatory  flow  was  studied  by  Nakato  et  al  (1977)  and 
MacDonald  (1977) .  All  of  their  measurements  show  that  the  turbulence 
intensity  as  well  as  the  sediment  concentrations  decay  exponentially 
away  from  the  bed 
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c(z) 


(2.9) 


Figure  2.6:  Distributions  of  vertical  velocity  fluctuations ,  v 

ran 

and  tine  averaged  sediment  concentration*  c  from  MacDonald  (1977) . 
Both  are  exponentially  distributed  with  very  similar  length  scales 
lv  and  The  measurements  were  taken  over  a  circular  swing  moving 

through  3till  water.  7  *  6s,  a  *  0.27m,  w  «  O.OlSm/s. 
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Figure  2.7;  Tin*  averaged  suspended  sediment  concentrations  c  from 
Nakato  at.  al.  (1977),  measured  over  natural  ripples  in  an  oscillating 
water  tunnel.  T  ■  1.8s,  a  ■  .075m,  d  -  .14nm. 
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and  the  vertical  length  scales  are  very  similar.  See  Figures  2.6 
and  2.7.  This  is  a  very  important  observation  because  it  shows 
that  sand  and  turbulent  water  spread  away  from  the  bed  by  very 
similar  mechanisms  irrespective  of  the  fact  that  sand  ha.g  a  finite 
settling  velocity,  which  the  turbulent  water  has  not.  This  i3 
possible  because  of  the  way  in  which  sand  grains  travel  in  vortex 
flow.  We  shall  study  this  in  one  of  the  following  sections. 

Nalcato  et  al's  measurements  show  one  more  interesting  thing 
about  the  turbulence  intensity  over  ripples.  They  show  that  the 
intensity  is  increased  by  the  presence  of  moving  sand,  all  other 
things  being  equal.  This  is  surprising  and  so  far  unexplained , 
but  it  is  consistant  with  the  fact  that  flat  beds  of  loose  sand 
offer  much  more  resistance  to  oscillatory  flow  than  flat  beds  of 
fixed  sand  grains  (see  Figure  2.8) ,  and  also  with  the  fact  that 
Bagnold  (1946)  measured  no  energy  dissipation  factors  larger  than 
0.24  over  large  sharp  crested  ripples  in  clear  water,  while 
Carstens  et  al  (1969)  and  Lofquist  (1980)  observed  values  in  excess 
of  0.30  in  sediment  laden  water. 


» 


QUANTITATIVE  DESCRIPTION  OF  BOUNDARY  LAYERS 

We  shall  now  try  to  give  a  mathematical  description  of  the  flow  * 

in  the  wave  boundary  layer  that  isr within  four  or  five  ripple  heights 
or  a  few  hundred  grain  diameters  from  the  bed. 


We  assume  that  the  velocity  outside  the  boundary  layer  is  uniform 


1 


0 


in  spaca  and  varies  as  a  simple  harmonic  with  time 

u  *  u(t)  *  aiu  cos  iut  (2.11) 

OB  OB 

In  order  to  simplify  the  mathematical  treatment,  it  is  convenient  to 
use  the  complex  formalism 

: 

u^t)  -  cun  e1U)t  (2.12) 

I 

The  complex  number  cun  e1UJt  moves  around  a  circle  with  radius  cun  and 
angular  velocity  <u  while  the  real  part  cun  cos  tut  which  corresponds 
to  the  physical  velocity  oscillates  between  atu  and  -cun  on  the  real 

i 

axis,  and  the  imaginary  part  i  cun  sin  at  oscillates  between  the 
corresponding  values  on  the  imaginary  axis,  see  Figure  2.9. 

The  equation  of  motion  is 


3u  m  1  3t  1  3p 

3t  ”  p  3z  ”  p  3x 


(2.13)  : 

I 

i 


where  p  is  the  pressure.  This  equation  is  valid  both  outside  and 
within  the  boundary  layer  but  outside  the  boundary  layer  there  are* 
I  by  definitioth.no  shear  stresses  in  a  purely  oscillatory  flow.  We 

have  therefore 


3uo, 

3t 


1  3n 
p  3x 


(2.14) 


This  expression  for  the  pressure  gradient  is  valid  inside  the  boundary 


layer  as  well  if  the  pressure  distribution  is  hydrostatic.  That  is , 


aojcoscjt 


iuit 

Figure  2.9:  The  complex  velocity  aue  has  constant  modulus 
can  and  neves  around  a  circle  with  angular  velocity  ui.  The  real 
part  can  cos  ait  which  represents  the  physical  velocity  os¬ 
cillates  between  atn  and  -can. 


if  vertical  accelerations  in  the  boundary  layer  are  small  compared 
to  gravity.  In  that  case  equation  (2.13)  can  be  written 


_3_ 

3t 


(u  -  u) 
00 


_l 

p  3z 


(2.15) 


It  is  convenient  now  to  write  the  velocity  deficit  u  -  u  in  the 

00 

form 


(u^-  u)  *  <201  D  (z,t)  (2.16) 

where  the  non  dimensional  D(z,t)  shown  in  Figure  2. 10* is  a  complex 
function  which  must  be  equal  to  eia,t:  at  the  boundary  where  the  velocity 
vanishes 

iujt 

u(0,t)  »  0  -*■  D(0,t)  -  e  (2.17) 

At  infinity  the  deficit  is  of  course  zero  so  we  also  have 

D(z,t)  -*■  0  for  z  (2.18) 

The  deficit  function  D(z,t)  takes  a  particularly  simple  form  when 
the  bed  is  smooth  and  the  flow  is  laminar.  For  laminar  flow  the 
shear  stress  is  proportional  to  the  velocity  gradient 

3u 

T  -  P  V  —  (2.19) 


where  vis  the  kinematic  viscosity  of  water,  so  equation  (2.15)  can 
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be  written 


_3_ 

3t 


[aiuD  (z  ,t)  ] 


1  3 


[u^- atiiD(z,t)  ]  ) 


which  reduces  to 


3D  3*0 

3t  *  V  3  z* 


(2.20) 


(2.21) 


We  solve  this  equation  by  separation  of  the  variables,  z  and  t. 

First  we  notice  that  the  velocity  deficit  D(z,t)  must  be  periodic 

in  time  like  u  (t)  so  we  assume  the  following  form 
00 


D(z,t) 


t  D  (z) 
l  a 


incut 

e 


(2.22) 


which  we  insert  into  (2.21)  and  get 


inu  D 

a 


(2.23) 


This  has  the  solution 


D  (z)  -  A  e/inw/V  z  +  B  ,-/ina)/V  z 

h  n  n 


(2.24) 


Now  since  must  decrease  with  increasing  z  we  must  have 


0  for  all  n. 


(2.25) 
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which  leaves  us  with 


,  9  _  .  .  incut 

u(z,t)  ■  u  (t)  •  an  i  D  (z)  e 

*  in 


(2.26) 


icut 

*  atue  -  aw 


?  -/imu/v  z  incut 

o  3  Q 


l  n 


(2.27) 


The  velocity  must  vanish  at  the  bed  so  we  have 


0  -  ® u  ?  B  einUt 

l  n. 


(2.28) 


which  gives 


a  «  i 

i 


(2.29) 


0  Cor  n  *  1 


(2.30) 


so  Cor  laminar  Clow  over  a  smooth  wall  we  have 


D(*,t)  -  D  (zie**  -  e'-^'e^ 

l 


(2.31) 


and 


u(z,t)  -  ocu  Cl-  •“’W  2]  eiwt 


(2. 32) 


Although  laminar  Clow  over  a  smooth  bed  is  a  rather  special  case, 


unlikely  to  occur  in  nature,  the  solution  (2.31)  Cor  this  case  shows 
most  oC  the  general  properties  oC  oscillatory  boundary  layers. 
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If  ve  write  (2.32)  in  the  form 


u(z,t)  »  [1  -  D(z)l  aweiUt 

l 


(2.33) 


we  can  see  how  o  (z)  determines  the  local  amplitude  and  phase  relative 
to  the  free  flow. . 

We  see  that  when  D  (z)  is  a  real,  negative  number,  the  velocity 
is  in  phase  with  u  (t)  and  the  amplitude  is  larger  than  a us.  That 
is,  there  are  elevations  where  the  velocities  get  larger  than  those 
of  the  free  flow. 

Remembering  the  form  of  □  , 


D  (z) 
l 


(1+i) z 


(2.34) 


or 


D  (z) 
l 


|”cos 


-  isin 


/u>/2v 


J 


(2.35) 


« 


we  see  that  this  occurs  for 


C 


«  pir  ,  p 


1,3,5... 


(2.36) 


Figure  2.11a  &  b  show  the  position  of  D  z)  and  of  l-D(z)  in  the 
complex  plane.  Figure  2.11e  shows  the  position  of  u(z,t)  in  the 
complex  plane  and  Figure  2. lid  shows  the  variation  of  the  velocity 
amplitude  a(z)  with  the  non  dimensional  elevation,  z/gj/2V. 

The  general  picture  of  the  velocity  variation  is  the  same  for 


Figure  2.11;  Velocity  variations  with  elevation  in  oscillatory  laminar 
flow  over  a  smooth  wall. 

a:  The  deficit  function  D(z,t)  moves  along  a  logarthmic  spiral  starting 

at  1  and  approaching  0  as  z  increases.  Numbers  on  the  curves  refer 
to  the  non  dimensional  elevation  /w/2 Vz. 
b:  Corresponding  variation  of  l-D^Cz)  which  is  the  ratio  between 

u(z,t)  and  uM(t)  zee  equation  (2.33). 
c:  In  the  simple  case  of  laminar  flow  over  a  smooth  wall  where  u(z,t) 

is  simple  harmonic  (u(z,t)  s  awd-D^z)  )eiajt)  ,  we  can  construct 
u(z,t)  geometrically  by  using  the  circle  from  figure  2.9  and  the 
spiral  of  figure  2.11b). 

d:  The  variation  of  the  velocity  amplitude  U(z)  *  au\ 1-Di  (z) |  with 

elevation . 
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turbulent  flow  with  small  roughness  to  semi  excursion  ratios 
(r/a  £  0.1)  as  shown  by  Jonsson  (1980).  Again  the  path  of  D^(z) 
in  the  complex  plane  is  a  spiral  which  starts  at  1  for  z  *  0 
and  approaches  zero  quickly  as  z  grows.  The  expression  (2.34) 
shows  that  for  laminar  flow  over  a  smooth  boundary ,  the  complex 
logarithm  of  D  (z)  is  a  particularly  simple  function  of  z: 

In  D  -  -(1  +  i)  /372V  z  (2.37) 

l 

and  that  the  angle  Arg  (D^)  and  the  logarithm  of  the  modulus  In  |d^ | 
are  identical: 


\  Arg  (D  )  »  In  Id  I  -  -/w/2v‘  z  (2.38) 

l  '  x 

Kalkanis'  (1964)  measurements  show  that  for  turbulent  flaw  aver 
5  a  smooth  wall  the  relation  is 

Arg  (D  )  »  1.2  In  Id  |  (2.39) 

l  i 

I 

and  both  increase  as  z° • 3  approximately.  See  Figure  2.12.  For  flow 
over  two  dimensional  roughness  the  pattern  is  different  again. 

B  Figure  2.13  shows  experimental  data  from  Kalkanis  (1964)  ,  where  the 

bed  roughness  was  made  up  by  half  cylinders.  Figure  2.14  shows  data 
from  Jonsson  and  Carlsen  (1976)  who  worked  with  considerably  larger 
•  Reynolds  numbers  than  Kalkanis,  and  with  a  different  type  of  two 

dimensional  roughness  elements. 

The  data  plotted  in  Figures  2.12  through  2.14  show  some  variations 


in  behaviour,  however  a  general  feature  is  that  In  D  behaves  very 


Figure  2 . 12 :  Variation  of  InD.  for  turbulent  flow  over  a  smooth 
wall .  Kalkanis  (1964),  Run  105,  T  »  1.4s,  a  *  .25m. 


Figure  2.13:  Variation  o£  lnD^  for  flow  over  two  dimensional 
roughness  elements.  KaDcanis  *(1964),  Run  224,  T  »  1.4s, 
a  «  ..15m,  Jonsson  (1980)  estimates  that  the  roughness  in  Kalkanis 
experiments  was  equal  to  the  radius  of  Kalkanis '  roughness  cy¬ 
linders  ,  for  this  experiment  :  r  -  .010m. 
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Figure  2.14:  Variation  of  lnD^  for  rough  turbulent  flow.  Measurements 
from  Jonsson  and  Carlson  (19767.  The  values  obtained  above  z/<5i  ■  5 
are  probably  unreliable.  The  scaling  length  5i  is  the  real  displacement 
thickness  defined  in  equations  (2.59)  through  (2.61). 
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similarly  and  that  both  Arg(D  )  and  ln|D^[-  can  be  determined  by  one  power 
function  throughout  the  whole  experiment al  range.  This  is  interest¬ 
ing  because  all  of  the  theoretical  (or  speculative)  models ,  eg. 

Kajiura  (1968)  and  Jonsson  (1980)  work  with  two  or  more  layers  where 
the  behaviour  of  D  is  expected  to  be  different. 

l 

He  shall  not  however, go  deeper  into  the  study  of  the  details 
of  u  (z,t)  at  the  moment.  Instead  we  will  consider  another 
measurable  feature  of  the  boundary  layer,  namely  the  bed  shear 
stress  t  (0,t) ,  and  show  how  the  empirical  information  on  T  can 
be  used  to  predict  vertical  length  scales  in  the  boundary  layer. 

This  is  important  because  we  need  the  vertical  length  scales  for 
prediction  of  suspended  sediment  distributions  and  the  only  detailed 
empirical  information  about  the  boundary  layer  over  natural  sand 
bads  (Cars  tens  at  aL  1969  and  Lofquist  1980)  are  obtained  by  friction 
or  energy  dissipation  measurements ,  without  detailed  velocity  measure¬ 
ments. 

He  start  again  with  the  equation  of  motion 

3  ,  .  1  3t 

3t  '*•’  u  "  p  3z  (2.15) 

which  we  integrate  with  respect  to  z  and  get 

3  °® 

3t  {  ('W  u  >  dz  "  £  CT(0,t>  -  T(°°,t)  ]  (2.40) 

The  term  T(®,t)  can  be  neglected  since  the  shear  stresses  vanish 
outside  the  boundary  layer,  so  we  have 
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j 


t 


It  {  (U-'  U  }  dZ 


—  T(  0,t) 

P 


(2.41) 


which  shows  that  the  instantaneous  bed  shear  stress  is  proportional 
to  the  total  acceleration  deficit. 

Since  both  the  velocity  deficit  ua-  u  and  the  bed  shear  3tress 
T(0,t)  must  be  periodic  with  period  T  we  may  write  both  as  a  Fourier 
expansion 


u -  u 


aw  £  D  (z)e 
1  n 


inwt 


(2.42) 


t(0,t) 


p  (aai)  2  Z  C  e 
l  a 


inuit 


(2.43) 


After  introduction  of  these  and  cancellation  of  common  factors,  the 
integrated  equation  of  motion,  (2.41) ,  becomes 


intu  f  D  (z)dz 


aw  C 


(2.44) 


We  are  particularly  interested  in  the  fundamental  mode  (n  ■  1)  for 
which  we  get 


i  /  D  (z)  dz  «  a  C 
o  l  i 


(2.45) 


The  integral 


i  /  D  (z)  dz 

9  1 


(2.46) 
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is  the  complex  displacement  thickness  of  the  fundamental  mode  of  the 
velocity,  so  equation  6.45)  enables  us  to  determine  a  vertical 
length  3cale  for  the  boundary  layer  if  we  know  the  bed  3hear  stress 
and  vice  versa. 

We  shall  illustrate  this  by  looking  once  mere  at  the  special 
case  of  smooth  laminar  flow  where  we  have  found 


□  (z) 

l 


for  which 


A  *  i  f  D  (z)  dz  *  yiv/w 
l  •  l 


(2.34) 


(2.47) 


and  chus  from  (2.45): 


C  -  /iv/a^w  (2.48) 

l 


so  the  bed  shear  stress  in  smooth  laminar  flow  is  given  by 


T  (o,t) 


P  {cun) 


iajt 


(2.49) 


The  fact  that  the  argument  of  is  tt/4  means  that  the  bed  sheer 

stress  leads  the  free  flow  u#  «  m  exp  [iuit]  by  tt/4  or  4S  degrees. 

The  result  (2.49)  could  of  course  also  have  been  obtained  from  the 

definition  (2.19)  together  with  (2.34). 

Several  authors,  starting  with  Jonsson  (1966)  describe  the  bed 

shear  stress  in  terms  of  the  wave  friction  factor  f  defined  by 

w 
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ip  f  (aw) 3 

w 


(2.50) 


Normally  f  will  include  contributions  from  the  higher  harmonics 
(equation  2.43) 


■{;  •‘-J 


(2.51) 


and  there  is  no  3imple  relation  between  £  and  C  .  However  for 

w  x 

smooth  laminar  flow  where  only  C  is  non-zero  we  have 


if  -  C 
w  1  t 


(2.52) 


and  thus 


2/j/a2w 


(2.53) 


While  (2.52)  only  holds  for  smooth  laminar  flow  we  can  find  a  similar 


relation  for  the  energy  dissipation  coefficient:  f@,  which  holds 


as  long  as  the  free  flow  is  simple  harmonic  (u  (t) 


aw  exp  [ iwt ] ) . 


The  average  rats  of  energy  dissipation  per  unit  area  over  one 


period  is 


u„(t)  T(0,t) 


(2.54) 


see  eg.  Kajiura  (1968).  Since  ua(t)  is  assumed  simple  harmonic,  only 
the  fundamental  mode  of  the  shear  stress  will  contribute  to  the 


average 
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Re  |  aw  eiajt  }  Re  £  p  (aw)  *  C  eiUt  J 


E  -  i  p(aw) 3  Re  (C  ) 

l 


(2.55) 


(2.56) 


The  traditional  way  of  describing  the  energy  dissipation  is  in  terms 
of  the  energy  dissipation  factor  ,  defined  by 


E  -  f9  (au>* 


(2.57) 


which  is  related  to  C  by 


X 


f.  ■  T  R*{c1} 


(2.58) 


¥ 


As  mentioned  above,  most  of  the  existing  empirical  data  on  oscillatory 

boundary  layers  over  natural  sand  beds  consist  of  friction  or  energy 

dissipation  measurements ,  so  we  can  predict  f  or  Re  f  C  1  reasonably 

e  *■  l J 

well  using  those  data. 

In  order  to  utilize  the  same  data  for  the  prediction  of  vertical 
length  scales  in  the  boundary  layer,  we  now  look  for  one  which  is 
directly  related  to  or  Re  fcJ  ■  Equation  (2.45)  shows  that  what 
we  are  looking  for  is 


Re 


■£i  /  (z)  dz  J 


(2.59) 


which  is  related  to  by 


(2.60) 
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and  to  f  via  (2.58) 
s 


6  -  f  <2  (2.61) 

1  3tt  9 

PREDICTION  OF  f  ,  f  and  6, 
e  w  i 

Both  the  wave  friction  factor  and  the  energy  dissipation  factor 
will  in  general  depend  on  the  flow  Reynolds  number,  a2 oj/v,  as  well 
as  on  the  roughness  to  semi  excursion  ration,  r/a?  but  when  the 
flow  is  rough  turbulent ,  as  is  the  case  under  field  conditions  we 
can  reglect  the  Reynolds  number  dependence.  The  hydrodynamic  conditions 
of  roost  wave  flume  experiments  are  in  the  region  of  transition  between 
laminar  and  turbulent  flow,  see  Jonsson  (1966).  For  such  conditions 
there  is  no  empirical  data  (known  to  the  author)  about  friction  factors 
over  natural  sand  beds.  The  following  derivations  are  therefore  re¬ 
stricted  to  rough  turbulent  flow  conditions. 

The  wave  friction  factor  f^  and  the  energy  dissipation  factor 
fQ  are  theoretically  different  and  may  indeed  take  very  different  values 
under  some  conditions  (Figure  2.15).  The  differences  are  however, 
usually  small  compared  to  the  scatter  in  the  empirical  data.  There¬ 
fore  it  is  not  really  worth  while  to  distinguish  between  fQ  and  f^ 
for  predictive  purposes. 

In  the  following  we  will  make  use  of  data  from  Kamphuis  (1975) 
who  measured  shear  stresses,  and  thus  f^,  over  fixed  sand  bads,  from 
Carstens  et.  al.  (1969)  who  measured  energy  dissipation  (»  ffl) 
over  natural  sand  beds,  and  from  Lofquist  (1980)  who  measured  both, 
over  natural  rippled  sand  beds.  The  differences  between  and  f^ 
are  illustrated  by  Lofquist 's  data  in  Figure  2.15.  Only  two  experiments 
show  a  significant  deviation  corresponding  to  a  particularly  peaked 
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» 


I 
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Figure  2.15:  Comparison  of  measured  values  of  and  fe 

from  Lofquist  (1980.  We  see  that  and  fQ 
are  generally  nearly  equal.  The  few  extreme 
deviations  stem  from  experiments  where  r(t) 
has  a  pronounced, narrow  peak. 

variation  of  T(t).  The  particular  experimental  conditions  that  lead 
to  such  behaviour  are  as  yet  unknown ,  but  they  may  become  clear  with 
the  progression  of  Lofquist's  work. 

Prediction  of  ffl  for  rough  turbulent  oscillatory  flows  involves 
two  steps,  namely,  estimations  of  the  hydraulic  roughness,  r,  of  the 
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bed  and  determination  of  f  from 

e 


•exp  [5.213  (-)*1*  -  5.977] 

a 


for  | 


.63 


(2.62) 


0.30 


for  —  >  .63 

a,  — 


Equation  (2.62)  which  corresponds  to  the  flow  model ,  suggested  by 
Jonsson  and  Carlsen  (1976)  is  only  one  of  several  semi  empirical 
formulae  of  the  upper  limit  of  0.3  is  somewhat  arbitrary.  It  is 
also  doubtful  whether  all  of  the  underlying  physical  assumptions , 
such  as  the  existence  of  an  overlap  layer  (a  layer  where  the  shear 
stress  is  proportional  to  the  distance  from  the  bed  and  the  velocity 
gradient)  are  met  by  the  flow  over  large  natural  sand  ripples.  Me 
will,  however,  use  equation  (2.62)  indiscriminately  with  all  the 
experimental  data  from  Cars tens  et  al  (1969)  and  Lofquist  (1980) 


and  not  claim  that  title  resulting  predictive  formulae  are  anything 
but  empirical. 

The  bed  roughness,  r,  depends  on  the  bed  form  geometry,  sand 
grain  size  and  on  the  motion  of  sandgrains  near  the  sand-water- 
interface.  The  simplest  case  is  that  of  the  flat  bed  with  no  sedi¬ 
ment  motion.  Kamphuis  (1975)  measured  friction  factors  under  such 
conditions .  Figure  2.16  shows  that  his  results  are  well  predicted 
by  equation  2.62)  with 
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r 


2.5d 


(2.63) 
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Figure  2.16:  Friction  factors  measured  over  flat  beds  of  fixed 
and  loose  sand.  The  moving  sand  over  a  loose  bed 
will  increase  the  roughness  by  one  or  two  orders 
of  magnitude,  depending  on  the  intensity  of  the 
sediment  motion. 

When  the  sand  is  allowed  to  move,  the  friction  is  increased 
very  considerably,  3ee  Figure  2.16.  Grant  and  Madsen  (1982)  suggested 
that  the  roughness  contribution  from  moving  sand  be  given  by 

§■  -  160  9  (/I7  -  0.7  /F)1  (2.64; 

a  c  c 

which  is  based  on  the  mechanics  of  sand  saltation  in  air.  It  is 
reasonable  to  expect  r  to  be  a  function  of  the  shin  friction  Shields 
parameter  9'.  However,  the  form  of  (2.64)  seems  not  to  work  too  well 
for  the  data  of  Carstens  et.  al.  (1969) ,  see  Figure  2.17.  The 
expression  (2.64)  varies  essentially  as  9 ' f where  as  the  data  indicate 
a  slower  growth,  more  like  /S"5*  for  9*  >  0.1. 


Figure  2.17:  Observed  roughness  (via  equation  2.62)  for 

flat  beds  of  loose  sand.  Measurements  by 
Carstens  et.  al.  (1969) . 

When  the  bed  is  covered  by  ripples,  these  will  contribute  to 
the  roughness  as  well.  Grant  and  Madsen  argue  that  as  a  first 
approximation,  this  contribution  should  be  proportional  to  the 
ripple  height  n  and  the  steepness  n/X 
The  equation 

r  ■  9  n2/X  +  190  /9'  -  0.05'  d  (2 

provides  a  good  fit  to  all  the  experimental  data  from  Carstens  et. 
al.  (1969)  and  Lofquist  (1980).  See  Figure  2.18. 

The  skin  friction  Shields  parameter  9’  expresses  the  ratio 
between  the  moving  shear  force  (*•  Td2)  and  the  stabilizing  gravity 
force1  (-  pg  (s  -  1)  d3).  Zt  may  be  estimated  by 
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Figure  2.18:  Predicted  (eq.  2.65)  and  observed  (via  eq.  2.62) 
values  of  r/a  for  the  energy  dissipation  measurements  of  Carstens 
et.  al.  (1969)  A:  rippled  bed,  d  »  .19sn.  8:  rippled  bed,  d  - 
.30mm.  C:  rippled  bed,  d  *  .59m.  F:  flat  bed,  d  *  .30m  and  from 
Lofquist  (1980),  L:  rippled  bed,  d  -  .55mm.  The  straight  line 
corresponds  to  perfect  agreement  between  prediction  and  observation. 
The  corresponding  f9  values  will  be  in  closer  agreement  because 
varies  slowly  with  r/a. 


where  f'  is  calculated  from  (2.62)  using  the  roughness  of  a  flat 
w 

bed  of  fixed  grains:  r  •  2.5d,  s  is  the  relative  density  of 

the  sand  and  g  is  the  acceleration  of  gravity. 

For  prediction  of  the  ripple  geometry  under  field  conditions, 
one  may  use 


n/X  »  0.342  -  0.34  VF 


for  9*  <  1 


(2.67) 
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-um 


aa  suggested  by  Nielsen  (1981).  When  0’  >  1,  the  ripple  height 

should  be  taken  as  zero. 


TURBULENCE  STRUCTURE  OP  WAVE  BOUNDARY  LAYERS 

Oie  turbulence  strusture  of  oscillatory  boundary  layers  can  1 

be  observed  in  different  ways.  Indirectly  by  determinations  of  the 
eddy  viscosity  from  velocity  measurements 

» 


v 


T 


h  i  (u  ~  ucc}  d* 

3u  ~ 
3z 


(2.69) 


» 


as  done  by  Horikawa  and  WatanaJbe  (1968) ,  Lundgren  (1972) ,  Jonsson 

and  Caxlsen  (1976) ,  or  by  direct  measurement  of  the  vertical  velocity  I 

fluctuations  (v  )  as  done  by  Nakato  et.  al.  (1977) ,  MacDonald  (1977) , 
as 

Kemp  and  Simons  (1982)  and  Du  Toit  and  Sleath  (1981) . 

Local  eddy  viscosities  determined  by  (2.69)  taking  T  (z,t)  and  p 

u  (z,t)  as  the  real  physical  quantities  show  a  strong  variation 
with  tiem,  which  is  at  first  sight  quits  astonishing.  It  behaves 
like  -tan  (wt) ,  and  tends  to  minus  and  plus  indinity  on  either  p 

side  of  two  vertical  asymptotes  during  each  wave  cycle,  see  Figure 
2.19.  The  asymptotes  are  of  coarse  at  the  phases  where  the 

velocity  gradient  is  zero.  Now  the  eddy  viscosity  is  normally  » 


interpretatad  as  the  product  of  a  typical  turbul uent  velocity  and  a 
typical  length 


(2.70) 


It  is  not  physically  meaningful  however,  to  interpret  the  infinite 
or  negative  values  of  v  as  a  result  of  the  negative  velocity  - 
or  length  scales,  the  alternative  is  to  interpret  as  a  complex 
number  which  has  a  constant  modulus 


(2.71) 


and  an  argument  which  is  equal  to  the  phase  difference  between 


Figure  2.19;  Eddy  viscosities  from  Horikawa  and  Watanabe  (1968) , 
tained  as  instantaneous  ratios  between  shear  stress  and  velocity 
gradients.  The  vertical  assympthotes  are  at  the  phases  where  the 
velocity  gradient  is  zero. 
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shear  stess  and  velocity  gradient.  If  both  the  shear  stress 
T  (2ft)  and  the  velocity  gradient  vary  as  simple  harmonics 
with  time  we  may  describe  them  by  the  real  parts  of  the  following 
complex  functions. 


3u 
3  2 


T 


0‘  (2) 


t(2)  a 


iait 


e 


i(ut  +  $) 


(2.72) 


(2.73 


where  D  is  the  local  velocity  amplitude  and  the  prime  denotes 
differentiation  with  respect  to  z  and  t(z)  is  the  local  <*h e»r 
stess  amplitude.  The  corresponding  complex  eddy  viscosity  is 


v  -  ¥2.  -  Hal- .  ei4> 

t  3u  pa  •  (2) 

32 

1 


(2.74) 


which  is  a  complex  constant  (for  fixed  z)  with  modulus 


i 


T(Z) 

P  0* (2) 


(2.75) 


and  argument  $. 

In  contrast ,  an  apparently  time  dependant  eddy  viscosity  results 
if  we  take  the  ratio  of  the  real  parts  of  (2.72)  and  (2.73): 


Re/ 1  (z)  e1^***) 
pRe  {a*  (z)  elut  } 


(2.76) 


I 


» 


T  (z)  COS  (Olt  +  (&) 
0  u '  (z)  cos  ut 


(2.77) 


1 

4 


T  (z) 

0  U'  (Z) 


[cos  <j>  -  sin  $  tan  wt] 


(2.78) 


which  describes  the  curves  shown  in  Figure  2.13. 

The  fact  that  the  apparent  eddy  viscosities  measured  by  Horikawa 
and  watanabe  approach  the  asymptotes  from  the  left  through  negative 
values  means  that  sin  <J>  is  positive  (3ee  equation  2.78)  and  thus  that 
the  argument  <J»  of  the  complex  eddy  viscosity  is  positive.  For  the 
experiments  of  Jonsson  and  Carlsen  (1976) ,  $  is  always  less  about 

5  degrees  and  mostly  positive,  but  for  flow  over  fully  developed  sand 
ripples  it  might  be  larger  and  more  variable  due  to  the  strong 
vortices  which  may  cause  a  substantial  momentum  transfer  (*  shear 
stress)  even  when  the  velocity  gradient  is  zero. 

Lundgren  (1972)  derived  the  magnitude  |v^(z)  [  from  the  measure¬ 
ments  subsequently  published  by  Jonsson  and  Carlsen  (1976)  and  re¬ 
commended  the  following  formula  for  prediction  of  |vT|. 


1  t  1.34  wShT  £  exp[£] 


(2.79) 

'  1 


where  the  boundary  layer  thickness,  6,  is  defined  as  the  lowest 
elevation  where  the  velocity  amplitude  equals  that  of  the  free  flow 
(C(5)  •  am) .  Following  (2.79)  |v^|  will  increase  linearly  with  z 

for  z  «  5  and  decrease  exponentially  with  z  for  z  »  5.  Figure  2.20 
shows  a  somewhat  different  expression 
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0.4  u#  2 

V  -  -  (2.80) 

1  +  0.2  f^exp 

which  uses  the  displacement  length  51  for  scaling,  because  it  is 

predictable  on  the  basis  of  Carstens  et.  al,'s  and  Lofquist's  energy 

dissipation  measurements  over  natural  sand  beds. 

Figure  2.20  also  shows  measured  vertical  velocity  fluctuations , 

from  Nakato  et.  al.  and  MacDonald.  We  see  that  v  decays  ex¬ 
tras 

ponentially  like  jv^|  for  :  2  S,  but  at  a  somewhat  slower  rate: 

v  *  0.9  u#  a“z/4*5<5l  for  (2.81) 

ns  1 

For  z  <t  5.  ,  v  is  nearly  constant  at  least  half  way  down  to  the 
1  ms 

bottom  of  the  ripple  trough  (2  *  -n) . 


vrma  ~  °-65  u*  from  -n/2  <  2  S.  5l  (2.82) 


If  the  involved  experiments  are  comparable  in  spite  of  their  very 
different  Reynolds  numbers  and  r/a  values  we  may  obtain  an  estimate 
of  the  length  scale  l v  by  dividing  the  above  expressions  for  |v^| 


V 


(2.83) 


The  result,  shown  in  Figure  2. 2o,  should  not  however* be  taken  as  more 

than  an  indication  of  the  magnitude  and  variation  of  Experiments 

where  both  T,  r—  and  v _  are  measured  are  clearly  needed  for  a  deeper 

o  z  ms 

understanding  of  this  phenomenon. 
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Figure  2.10:  Turbulence  structure  of  oscillatory  boundary  layers. 

Left:  Vertical  velocity  fluctuations  measured  by  Macdonald  (•),  over 
artificial  roughness ,  made  up  by  half  cylinders  with  diameter  .03m; 
a  »  .282m,  T  -  6s,  r  -  0.015m,  and  by  Nakato  et.  al.  (•■* :  over  ripple 
crest,  at  over  -trough)  over  natural  ripples,  a  m  .075m,  T  ■  1.8s, 
d  *  0.14mm. 

Center:  Inferred  turbulence  length  scale  (eq.  2.83) . 

Right:  |v  j  from  measurements  by  Jonsson  and  Carlsen,  derived  by 

harmonic  analysis  and  equation  (2.75).  TEST  I:  a  *  2.85m  ,  T  *  8.39s 
r  -  0.023m.  TEST  II:  a  -  1.79m,  T  ■  7.20s,  r  -  0.062m. 
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WAVE— CURRENT  BOUNDARY  LAYERS. 

The  water  motion  in  the  surf  zone  normally  involves  oscillatory 

components  of  very  different  periods ,  from  the  order  of  one  second 

for  the  harmonics  and  incoming  waves  to  several  hours  for  currents 

that  are  forced  by  the  tides.  The  ability  of  these  components  to 

penetrate  to  the  bottom  and  contribute  to  the  motion  of  the  large  ' 

near  bed  sediment  concentrations  is  highly  variable.  The  short 

period  oscillations  develop  thin  boundary  layers  and  the  long  period 

oscillations  develop  thick  boundary  layers.  The  thickness  5  of  an  ► 

oscillatory  boundary  layer  is  proportional  to  the  square  root  of  the 
period: 

5  ^  /T  ~VT'  (2.84)  ) 

where  T  is  the  period  and  is  the  typical  eddy  viscosity.  This 

means  that  if  the  boundary  layer  related  to  a  ten  second  swell  wave  > 

is  S  centimetres  thick,  then  the  surf  beat  current  with  period  90 

seconds  will  have  a  boundary  layer  15  centimetres  thick  and  a  tidal 

current  with  period  12  hrs  will  have  a  boundary  layer  thickness  of  I 

about.  10  metres.  So  under  normal  surf  zone  conditions  with  depth  of 

the  order  2  metres P  the  wave  and  surf  beat  boundary  layers  are  thin 

conpared  to  the  depth  while  the  tidal  current  has  a  boundary  layer  9 

thickness  much  larger  than  the  water  depth. 

The  Figures  2.21  and  2.22  show  two  examples  of  surf  zone 
currents.  Both  were  measured  in  the  trough  behind  a  well  developed  ^ 

bar  where  the  depth  over  the  bar  crest  was  about  half  a  metre  and 


Figure  2.2  :  shore-normal  current  time  series  from  4  elevations  in  the  bar 
trough  showing  vertical  segregation: onshore  velocities  near  the  surface  and 
offshore  closer  to  the  bed.  An  electronic  filter  with  a  100  second  time  con¬ 
stant  has  been  applied  to  remove  incident  waves,  the  water  depth  was  between 
1.5  and  2  metres. 


LONGSHORE  VELOCITY  PROFILE  ,  MAY.  1982 
EASTERN  BEACH. LAKES  ENTRANCE, GIPPSLAND. 


A  10  SEC  FILTER. SAMPLED  EVERY  6  SECS. 
»  100  SEC  FILTER. SAMPLE  INTERVAL  1  MIN. 
□  100  SEC  FILTER. SAMPLE  INTERVAL  1  MIN. 
X  FAST.SAMPLE  INTERVAL  6  SECS. 


Figure  2.22;  Vertical  profile  of  tine-averaged  longshore  current 
velocity,  as  recorded  within  the  trough  of  Eastern  Beach  in  a  suc¬ 
cession  of  runs.  The  position  of  the  lowest  flow  meter  relative  to 
the  ripples  is  shown,  h  »  1.14m,  T  »  3.65s,  H  -  0.39m, 

n  -  0.08m. 


I 
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the  maximum  depth  in  the  trough  about  2  metres.  For  both  the  shore 
parallel  current  and  the  shore  normal  one,  we  see  that  the  mean 
(obtained  by  averaging  over  about  20  minutes)  varies  considerably  over 
the  depth.  The  shore  normal  mean  current  is  segregated  with  onshore 
motion  at  the  top  and  offshore  motion  near  the  bed;  while  the  long* 
shore  current  shows  the  usual  monotonous  increase  towards  the  sur¬ 
face.  the  wave  and  surf  beat  oscillations  however, are  unchanged 
right  down  to  the  lowest  flow  meter. 

The  shape  of  the  longshore  current  profiles  in  Figure  2.22  are 
remarkable  in  that  the  velocity  gradients  between  7  and  22  centimetres 
are  much  smaller  than  those  between  22  and  100  centimetres.  This  is 
because  the  waves  cause  a  vast  increase  in  the  eddy'  viscosities  in 
the  lower  layer.  These  eddy  viscosities  are  visualized  by  the  en¬ 
trained  sand  in  the  lee  vortices  shed  by  the  ripples.  The  particular 
measurements  shown  in  figure  2.22  have  been  discussed  in  great  detail 
by  Wright  et.  al.  (1982) . 

Detailed  mean  current  profiles  from  under  a  combination  of 
waves  and  currents  are  rare  and  the  laboratory  measurements ,  e.g. 
Brevik  and  Aas,  (1978)  and  Kemp  and  Simons  (1982)  are  all  performed 
with  the  current  parallel  to  the  wave  propagation.  However,  all  the 
laboratory  measurements  and  the  field  measurements  of  Cacchioneand 
Drake  (1982)  show  that  the  steady  current  profile  is  logarithmic 
above  the  wave  boundary  layer  while  the  velocity  gradients  inside 
the  boundary  layer  are • suppressed  by  the  extra  turbulence  created 
by  the  waves.  The  resulting  appearance  of  the  current  profile  is 
illustrated  in  Figure  2.23  where  we  have  assumed  a  particularly 
simple  eddy  viscosity  distribution ,  which  was  suggested  by 


(c)  MEAN  LONGSHORE  CURRENT  (d)  MEAN  LONGSHORE  CURRENT 


Figure  2.23s  Characteristics  of  the  mean  current  profile  in  the  pre¬ 
sence  of  waves,  strong  enough  to  dominate  the  near  bed  turbulence  struc¬ 
ture. 

a:  The  time  averaged  shear  stress  is  linearly  distributed  and  may  be 

assumed  constant  for  *  <  l  if  l  «  h. 
b:  The  eddy  viscosity  distribution  is  parabolic  for  z  >  l  and  is 

a  constant  for  z  <  l. 

ci  The  u(z)  distribution  is  logarithmic  for  z  >  l.  With  no  waves  but 
the  same  u#,  the  profile  would  have  the  same  shape  (same  3u/3z) 
but  be  shifted  to  the  right. 

d:  Same  u(z)  profile  in  the  usual  semi  lagarithmic  plot  where  the 

logarithmic  part  becomes  a  straight  line. 
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Christoff ersen  1980. 


J 


v_  =  v  for  z  <  l 

T  w 


VT  =  tcu^z  (1-z/h)  for  z  Z  1 


(2.85) 


I 


I 

A 


The  parabolic  part,  used  for  :  >  I  is  the  usual  egression 
for  a  steady  flow  eddy  viscosity,  <  being  von  Kansan's  constant 
(«  0.4)  and  u#  being  the  time  averaged  friction  velocity  defined 
by 


u#  -  /T~(0)/f  1 


(2.86) 


where  x  (z)  is  the  time  averaged  shear  stress.  The  distribution  of 
X(z)  must  be  linear  and  we  neglect  wind  shear  stress  (x  (h)  a  0)  so 
we  have 


I 


x  (z)  -  x(0)  (1  -  z/h) 


or  in  terms  of  the  friction  velocity 


(2.87) 


X(z)  -  pu*  (1  -  z/h)  (2.88) 

When  the  distributions  of  shear  stresses  and  eddy  viscosities  are 
known,  we  get  the  velocity  gradients  3u/3*  from  the  definition  of 
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(2.89) 


With  the  expressions  (2.85)  and  (2.88)  this  yields 


3u 

3z 


u*  (1  -  z/h) 


v 

w 


for  z  <_  l 


(2.90) 


aS  , 

37  <  z 


for  z  >,  l 


(2.91) 


we  can  find  u(z)  for  z  <.  I  by  integrating  (2.90)  and  using 
u(0)  -  0: 

u2 

•  **#  "I 

u(z)  -  v""  (*  ~  2h  z2)  for  2  £  1  (2.92) 

w 

the  velocity  above  z  -  t  is  found  by  integrating  (2.91) 

«*  - 

u(z)  -  —  In  (-H  for  z  >  l  (2.93) 

K 

zx  represents  the  constant  of  integration  and  is  equal  to  the  level 
where  the  extrapolation  of  (2.93)  intercepts  the  z  axis,  see  Figure 
2.23c  and  2.23d.  The  value  of  zx  can  be  determined  by  matching  the 

two  expressions  (2.92)  and  (2.93)  at  z  -  l.  We  find 

zl  •  l  expC - - —  (1  -  (2.94) 

w 


> 


I 


» 
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or  for  Z  «  h 

<u#l 

zx  -  l  exp  [ - - — ]  (2.95) 

w 

Cacchoine  and  Drake  (1982)  found  zi  -  4.8  an  under  field  conditions 
and  Kemp  and  Simons  (1982)  found  zx  values  between  1.3  and  1.9  milli¬ 
metres  under  1  second  waves  with  am  increasing  from  0.06  to  0 . 12m/s , 
and  a  depth  averaged  current  of  0.184m/s.  The  bed  roughness  elements 
were  triangular  wooden  strips  of  height  5  millimetres,  and  the  z 
intercept  z„,  for  steady  current  alone  was  0.84  millimetres. 

Experimental  determination  of  \)^  is  difficult  but  Kemp  and 
Simons  found  that  the  assumption  of  being  a  constant  in  the  lower 
layer  was  quite  well  justified  under  their  .  tperimental  conditions 
where  am  was  of  the  order  0.3  to  0.7  of  the  depth  averaged  mean 
current. 

Brevik  and  Aas  (1980)  performed  similar  experiments,  only  with 
a  stronger  mean  current  relative  to  am  and  with  sinusoidal  ripples  1 

that  would  have  been  less  efficient  vortex  makers  than  the  triangles 

used  by  Kemp  and  Simons.  Their  velocity  profiles  were  logarithmic  ,j 

in  the  inner  layer  which  indicates  a  linearly  increasing  as  assumed 
by  Grant  and  Madsen  (1979) . 

The  limited  experimental  evidence  seems  thus  to  indicate  that  the 
ripple  eddy  viscosity  model  (2.85)  applies  when  the  boundary  layer  is 
dominated  by  the  waves, while  a  different  model  like  the  one  of  Sraat 
and  Madsen  (1979)  must  be  applied  when  the  current  dominates.  Visual 
observation  of  bed.  forms  can  easily  reveal  which  one  is  dominant 


at  a  given  location.  If  the  bed  is  flat  or  covered  by  symmetrical 
ripples,  it  is  nest  likely  that  the  waves  are  dominating;  but  if 
the  bed  is  covered  by  asymmetrical  dunes  like  in  rip  currents,  the 
current  must  be  dominating. 
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3 . _ THE  MOTION  OF  SUSPENDED  PARTICLES 


Since  the  main  topic  of  the  present  report  is  suspended  sediment  we 
shall  study  the  small  scale  motion  of  such  particles  in  some  detail. 

We  do  that  in  terms  of  an  analytical  solution  for  the  equation 
of  motion  of  a  suspended  particle  in  an  accelerated  non  uniform  flow 
field.  This  solution  treats  the  effects  of  gravity  and  those  of  the 
fluid  accelerations  separately,  which  is  a  major  advantage  because 
it  enables  a  step  wise  analytical  treatment  of  the  problem  which  is 
easier  to  follow  physically  than  numerical  models  that  take  all 
steps  at  once. 

Wave  boundary  layers  are  dominated  by  vortices  and  therefore 
we  pay  special  attention,  to  this  kind  of  flow  structure.  However 
the  resulting  formulae  are  quite  general  and  apply  to  any  flow  as 
long  as  the  fluid  accelerations  are  reasonably  small  compared  to 
the  acceleration  of  gravity. 

SEDIMENT  MOTION  IN  VORTEX  FLOW,  QUALITATIVELY 

For  kinematic  as  well  as  for  dynamic  analysis  of  suspended 
sediment  motion,  it  is  convenient  to  split  the  sediment  particle 
velocity  u^  into  three  parts,  as  shown  in  Figure  3.1 


Figure  3.1:  The  sediment  velocity  is  split  into  three  parts: 

-*■  -*■ 
u  ■  u  +>  w  +•  v,  and  v  is  treated  as  a  perturbation  to  u  +•  w. 

P 


u  (x,z,t)  «  u(x,z,t)  +  w.  ♦  v(x,z,t)  C 

P 

where  u  is  the  water  velocity,  w  is  the  still  water  settling 
velocity  and  v  is  due  to  the  accelerations  of  the  flow  (see 
Figure  3.1).  In  most  practical  situations  the  fluid  accelerations 
are  much,  smaller  than  the  acceleration  of  gravity,  g,  so  we  have 


and  can  use 


u  *  u  +  w  ( 

p 


as  an  approximation  for  the  sediment  particle  velocity.  Let  us 


see  what  that  leads  to  in  vortex  flow. 
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Consider  a  forced  vortex  with  the  velocity  field 


Q("Z) 

x 


(3.4) 


which  is  that  of  a  rigid  body  rotating  with  angular  velocity  (2. 
Applying  equation  (3.3)  we  find 


a(«  :sw/a> 


(3.5) 


(w  *  (0 ,  -w) ) .  This  shows  that  a  sediment  particle  will  move  along 

some  circle  with  centre  (w/Qf0)  (see  Figure  3.2).  Tooby  et.  al. 
(1977)  showed  experimentally  that  small  particles  do  in  fact  follow 
these  circular  paths  very  closely  and  only  very  slowly  spiral  away 
from  them.  He  shall  study  this  spiralling  process  in  detail  later. 

The  interesting  implication  of  these  approximatly  closed 

paths  is  that  the  Lagrangian  mean  velocity  for  the  sediment  particles 

is  approximately  0  and  not  w  as  one  would  find  from  u  »  u  -t»  w, 

P 

♦  + 
assuming  that  u  is  random  turbulence  with  average  0.  This  asphasizes 

the  importance  of  considering  the  flow  structure  in  the  context  of 

sediment  entrainment. 

The  fact  that  the  settling  velocity  is  partly  eliminated  in 
a  vortex  flow  field  has  been  touched  upon  by  Seizes  (1977)  who 
found  it  in  a  numerical  study.  Be  concluded  that  the  sediment  par¬ 
ticle  must  spend  more  time  in  the  upward  moving  parts  of  the  flow 
than  in  the  downward  moving  parts.  This  is  indeed  the  case;  the 
sediment  path  sketched  in.  Figure .3.2  Lias  completed  in  the  "upward 


Sediment  oath 


Woter  particle  path 


Figure  3.2:  Sediment  particle  with,  settling  velocity  w  moving  in  a 
forced  vortex  with  velocity  field  u.  Thfc  particle  paths  corres¬ 
ponding  to  u  +  w  are  circles  and  thus  the  lagrangian  mean  of  the 
sediment  particle  velocity  will  be  zero. 


moving1*  part  of  the  vortex.  For  buoyant  particles  or  air  bubbles 
the  center  of  the  "sediment  path"  would  lie  on  the  negative  x-axis 
and  thus  they  would  spend  the  majority  of  the  time  within  downward 
moving  fluid. 

the  described  trapping  mechanism  will  work  for  all  sand  grains 
with  settling  velocity  smaller  than  the  maximum  velocity  in  the  vortex 

the  next  questions  to  be  asked  is  whether  the  trapping  is  a 
feature  of  the  rather  unnatural,  forced  vortex  only.  The  answer 
is  no  for  the  following  reason  :  In  general,  the  velocity  field 
of  a  two  dimensional  vortex  can  be  written 
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u(X)  -  GF(x2  +  z2)("Z) 

z  x 


(3.6) 


and  che  first  approximation  to  the  sediment  velocity  is  then 


u  +  w  *  ftF(x2  +  z2) (  Z)  + 


(3.7) 


For  the  components  u  and  u  of  this  we  have 

px  pz 


V*’  - 21  -  -Vx'*’ 


and 


(3.8) 


» 


u  (X,  -  z)  -  u  (x ,  z) 
pz  pz 


(3.9)  » 


this  symmetry  means  that  any  particle  path  that  crosses  the  x-axis 
twice  must  be  closed  since,  due  to  the  syssnetry,  a  particle  which 
has  travelled  along  the  curve  PXP2  must  travel  bade  to  Px  via  the 
image  of  PXP2  (see  Figure  3.3) .  Hence  closed  sediment  paths  are  a 
general  feature  of  vortex  flow. 

A  fair  model  of  many  natural  vortices  is  the  Hankine  vortex 
in  which,  the  velocity  field  is  given  by 


-*■  ,x. 

ttV 


Or 

1  +  (x/R)2  +  (Z/R)2 


(3.10) 


I 


I 


Figure  3.5 1  Particle  paths  corresponding  to  u  »  u  +  w  in  a  Rankine 

P 

vortex.  In  the  inner  region  where  |u|  is  proportional  to  the  dis¬ 
tance  from  the  origin,  we  get  curves  that  are  very  similar  to  the 
circles  in  the  forced  vortex  (Figure  3.2). 


QUANTITATIVE  APPROACH 


Consider  a  particle  with,  diameter  d  and  relative  density  s, 
moving  in  a  liquid  with  velocity  u(x,2,t)  under  the  action  of 
gravity.  The  equation  of  motion  is: 
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The  first  term  on  the  right-hand  side  represents  the  pressure 
gradients  in  the  flow.  The  second  expresses  the  influence  of 
the  hydrodynamic  mass  which  is  C^s  times  the  mass  of  the 
particle.  The  third  term  is  due  to  the  drag  force  on  the  par¬ 
ticle,  and  is  the  instantaneous  drag  coefficient  which  is 
a  function  of  the  instantaneous  Reynolds  number. 


(3.14) 


The  validity  of  using  this  instantaneous  drag  coefficient 
C*  was  verified  by  Ho  (1964) ,  who  used  it  for  numerical  cal¬ 
culations.  In  the  following  we  use  the  approximation 


where  CQ  is  the  drag  coefficient  corresponding  to  the 
terminate  settling  velocity  through  still  water.  The  exponent 
y  is  thus  defined  by 


In  the  laminar  case ,  where  Stoke 1 s  law  is  valid  we  get 
y  -  1  and  in  the  turbulent  areas  where  CQ  is  practically  in¬ 

dependent  of  SI  we  get  Y  *  0.  See  Figure  3.6. 
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I 

Figure  3,6:  Hie  variation  of  the  drag  coefficient  CQ  with 
the  grain-Reynold 1 3  number  :  wd/v.  After  Navntoft  (1968) 


C 


The  Basset  term,  which  yields  a  drag  effect  when  |u-u  | 

P 

changes  magnitude  rapidly,  has  been  neglected  because 

[up-u(  -  (w  +  v(  (3.17) 

is  essentially  constant  (s  w)  under  the  assumption 


||/9<<1  (3.18) 

which,  holds  for  most  practical  cases.  From  the  definition  of  CQ 

(s-l) 0g  £  d3  PCQ  w2  J  d2  (3.19) 


we  get 


D 

4d(s-l) 


(3.20) 


Now  we  introduce 


u  (x,z,t) 
P 


u(x,z,t)  +  w  +■  v(x,z,t) 


(3.21) 


into  the  equation  of  motion  with 


3u  -*■ 

at  +  u‘7u 


(3.21) 


dv  a  v  . 

+  +  (u+w+v)  -Vv 


(3.22) 


then  we  get 


—  +  w*Vu  +  v*Vu  » 


(3.23) 


.1  .  ,du  M  d 

C?  1Jdt  *  T  dt  <— *’ 


The  drag  term  is  simplified  by  equation  (3.20),  and  we 


introduce 


(3.24) 


So  we  get  the  following  simplified  equation  of  motion 
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\ 


c 


S 


< 


dv  1  _ 

d t  1+cj/3  1+C^/s 


1  •* 
vVu 


du  g  I — *i  1-Y 

dt  ~i-7  iw+v'  (W+V)  "  g 

wz  T 


(3.25) 


To  assess  the  relative  magnitude  of  the  terms  we  bring  the  equation 
on  dimensionless  fora  by  the  following  transformations 


a 

9 

u/w 

(3.26) 

-♦> 

(3.27) 

V 

SB 

V/w 

-4. 

w 

w/w 

(3.28) 

T 

S 

1 

flt 

(3.29) 

X 

a 

X/R 

(3.30) 

z 

m 

X/R 

(3.31) 

Where. 

J2  is  a 

typical  angular  velocity  of  the  flow  and  R  a 

characteristic 

length.  This  leads  to: 

dV 

w/JlR 

S.n£  .  w/flR 

dT 

1+V3 

M  7u  *  l+ys  7  7u 

dO 

^dT  * 

&  l55?i1'Y  ^  -  A 

(3.32) 

I 


and  in  terms  of  the  vector  components  with  the  drag  term  expanded 


C 


« 


w/flR  3ax  w/qr  r  3UX 

1+C^/s  3  2  1+S/3  L*  3x 


V 

z 


(1-Y)  V  V  * 

X  2 


] 


(3.33) 
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£tt 

w/fla  _ 2  w/Ar 

l+Cj^s  3Z  1+C^/s 


v 

z 


We  omit  tarns  of  the  order  of  magnitude  |  v  [ 2  and  rearrange 

r 


dV 

X 

dT 

+ 

22.  v  + 

W (l  X 

w/Qr 

1+V“ 

1 — 1 

4s 

V 

X 

3  a  1 

+  TT  vzJ  “ 

da 

_ £  . 

w/OR 

3a 

X 

(3.35) 

* 

“0t 

dT 

1+V9 

3z 

dvz 

dT 

w/OR  j 

i+v*  L 

'3D* 

IT 

m 

3°z  1 

vx  *  7T  vzJ 

- 

dU 

W/  R 

1+V9 

3a 

-a 

15"  + 

z 

IT 

(3.36) 

1 

In  most  practical  cases  w  will  be  a  few  centimetres  per  second  or 
less  and  R  will  be  a  centimeter  or  more,  so  the  second  term  of  the 
left  hand  side  will  dominate  completely  over  the  third,  which  may  thus 
be  neglected.  The  resulting  aquations  on  dimensional  form  are 


When  the  velocity  £ield  is  steady  and  uniform,  the  driving 


terms  disappear  and  the  solutions  are 


Vel  -  W  •*»  [-  ^ 


V'1  . 


(t0)  exp  j^-(2-Y)^  tt-t0)  1 


which  shows  that  deviations  from 


u  ■  u  +  w  c: 

P 

decay  with  a  time  scale  of  which  is  for  most  natural  sediments 

og 

«•  2 

less  than  10  seconds. 

This  corresponds  well  with  the  results  of  Jenkin's  (1973)  who  found 
numerically  that  a  sand  particle  released  in  still  water  will  only 
travel  a  distance  comparable  to  its  own  diameter,  before  it  has  ob¬ 
tained  its  terminal  settling  velocity. 

Note  that  for  a  sand  grain  moving  in  air  the  situation  is 

w 

very  different  because  the  time  scale  —  is  at  least  one  order  of 

og 

magnitude  bigger  in  that  case  and  so  are  the  typical  velocities. 


In  an  accelerated,  flow  field,  deviations  from  (3.3)  will  be 

very  small  compared  to  w  because  of  the  large  damping  factor  22. 

(see  equation  3.37  and  3.38).  So  equation  (3.3)  can  be  expected  to 

give  a  good  first  approximation  to  the  sand  motion  in  a  flow  field 

i  du  i 

with  moderate  accelerations ,  (  -r—|  <  <  g) . 


SIMPLE  EXAMPLE 


To  gain  a  basic  understanding  of  the  acceleration  effects, 
let  us  look,  at  the  very  simple  example  of  a  small  grain  (Y»l) 
moving  in  a  forced  vortex.  The  velocity  field  is  given  by 


-►  x 

u$ 


(3.4) 


and 


u  +•  w  “ 


n<«0V 


(.3.5) 


Equation  (3.S)  corresponds  to  motion  along  any  circle  which  centre 
<w/0,0)  and  angular  velocity  ft-  See  Figure  3.2. 

For  our  small  grain  with.  y»l  equations  (3.37  and  3.38)  yield: 


&  *  sa  ;  .  -c 

dt  w 


d  f  *f!*\  .  w  3  /  -flz\ 

dt  ^  Ox  y  1+Cjj/s  Jz  V  Ox  / 


(3.41) 


We  introduce 


/  x\  m  /r  cos  Ot  +  w/fl 
\  zj  VR  0*  4 


(3.42) 


which  corresponds  to  the  same  steady  circular  motion  as  (3.5),  and 


t 
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1 

s 

* 


i 


(3.43) 


The  £irst  tern  on  the  right  hand  side  is  a  centrifugal  acceleration 
which  results  in  a  rotating  velocity  vR: 


fl2R  w  1 


rcos  (Qt  -  Tan”^  — ) 
ag 


-l  uG 

sin  (Qt  -  Tan  — -) , 

ag 


(3.44) 


In  most  practical  cases,  we  will  have  wfi/og  <  <  1  and  thus  we  see 

that  v  is  nearly  proportion  to  Q2R /g  -  j~-|/g  and  that  v  is 
R  at  R 

^  -t  tjfJ 

parallel  to  R  except  for  the  small  angle  Tan  — ■ -  (see  Figure  3.7). 

ag 

The  second  term  on  the  right-hand  side  of  (3.4  3)  is  the 

Coriolis  acceleration  which  gives  a  steady  drift  v  towards  the  left 

s 

for  s  <1  and  to  the  right  for  s  <1  (see  Figure  3.7) . 


(3.45) 


« 


The  effect  of  v  is  to  make  the  grain  spiral  out  (for 
w  >  0)  slowly.  To  estimate  the  time  scale  for  this  process  we  put 


v 


R 


Si  w  (COS 
g  sin 


nt 

nt* 


(3.46) 


which  assumes 


tgura  3.7; 


The  density  difference  causes  a  centrifugal  effect 


that  will  make  heavy  (s  > 1)  particles  spiral  out  and  light 
(s  < 1)  particles  spiral  in.  The  Coriolis  effect  leads  to  a 
steady  drift  to  the  left  for  s  >  1,  to  the  right  for  s  <1.  Note 

that  the  size  of  v  and  v  are  vastly  exaggrated,  compared  to 

-*■  -t.  K  * 

u  and  w. 


wfl 

og 


<  <  1 


and  write  v  in  the  form 
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The  last,  steady  term  is  not  contributing  to  the  spiralling 
process  in  average  over  a  revolution  so  the  spiralling  can  be 
modelled  by 


dR 

dt 


WR 


(3.49) 


(u+w  is  perpendicular  to  R)  which  has  the  solution 


R(t)  -  R„  exp  [  (t-t„)] 


(3. SO) 


showing  that  the  tima  scale  of  the  spiralling  process  is 


t  m  -2-  (3.51) 

3  w a1 

This  time  scale  will  normally  be  several  seconds  in  which  time 
the  particle  is  travelling  with  the  vortex  with  effective  settling 
velocity  zero. 

The  solution  given  by  (3.44)  and  (3.45)  is  valid  for 
particles  that  are  so  small  that  Stokes  law  applies,  that  is 
Y  ■  1.  The  general  solution  is: 

cos  (fit  -  Tan*1 6)  -  (3.52) 

is-ij  g 

vg  «  sintfit  -  Tan“lCS/(2-Y)  ])  (3.53; 


where  0  ■  wfl/ag  «  1  . 
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SECOND  ORDER  EFFECTS 

Due  to  the  non  linearity  (for  y  >  0)  of  the  drag  term 
in  equation  (3.13),  a  sediment  particle  will  settle  slower  through 
an  accelerated  or  non  uniform  flow  field  than  through  still  water. 

This  effect  can  be  analyzed  by  considering  second  order  terms 

* 

(o (V2 ) )  in  solutions  of  equations  (3.33)  and  (3.34).  Nielsen 

(1979)  studied  this  phenomenon  for  a  homogenous,  vertically  oscillating 

flow,  (u  »  Rflcosflt) ,  and  found  that  the  reduction  in  settling 
z 

2  2  > 
velocity  is  of  the  order  of  magnitude  [itfi/g]  w  which  is  in  most 

practical  cases  very  small  compared  to  w.  This  analytical  solution 

was  shown  to  be  in  good  agreement  with  measurements  and  numerical 

calculations,  carried  out  by  Ho  (1964).  * 

I 

» 

I 


4. 


FIELD  DATA 


THE  SEDIMENT  SAMPLER 

Over  the  twenty  month  period  from  October  1980  to  July  1982, 
the  Coastal  Studies  Unit  has  collected  71  detailed  profiles  of  time 
averaged  suspended  sediment  concentrations  c(z),  complete  with  hydro 
dynamic  data  and  bed  sediment  samples. 

The  sediment  sampler  which  was  designed  by  Peter  Nielsen  and 
constructed  by  Graham  Lloyd  is  probably  the  most  simple  instrument 
one  can  think,  of  for  this  kind  of  task.  It  has  no  electrical  or 
mechanical ,  moving  parts  and  is  therefore  very  rugged  and  dependable 
The  sand-water-mixture  is  collected  by  suction  and  the  driving  pres¬ 
sure  difference  is  provided  by  connecting  the  sample  jars,  that  sit 
on  the  sea  bed,  to  a  coianon  air  outlet  above  the  water. 


F: 


The  sediment  sampler 


The  sampler  is  designed  to  sample  a  vertical  array  of  simultaneous 
time  averaged  sediment  concentrations,  c(z). 

The  mav-i  nmm  distance  between  the  top  and  bottom  intakes  is 
.55m,  and  the  seven  intakes  can  be  placed  with  different  intervals. 

The  interval  between  the  two  bottom  intakes ,  where  the  concentration 
gradients  are  largest  is  normally  4  centimetres  and  the  intervals 
between  the  higher  intakes  vary  from  10  to  15  centimetres.  The  in¬ 
terval  between,  the  bed  and  the  bottom  intake  is  adjusted  once  the 
sampler  has  been  placed  in  the  sampling  position. 

The  intake  nozzels  are  stainless  steel  pipes,  15  centimetres 
long  and  with  inner  and  outer  diameters  of  3  and  6  millimetres 
respectively.  The  intakes  are  connected  to  the  sample  jars  by 
transparent  plastic  hose  with  5  millimetre  inner  diameter. 

Glass  3 airs  of  2.2  litre  capacity  are  used  as  sample  jars .  This 
size  is  sufficient  for  obtaining  measurable  amounts  of  sand  (2  .05 
grammes)  at  concentrations  as  small  as  10“ 5  by  volume.  Thi3  capacity 
is  normally  only  required  outside  the  surf zone.  For  measurements  in 
areas  with  breaking  waves  and/or  strong  net  currents  it  is  recommended 
to  use  smaller  jars,  which,  will  be  easier  to  handle. 

Because  the  jars  must  be  placed  on  the  bed  and  remain  stable, 
while  still  full  of  air,  a  considerable  mooring  weight  is  needed. 

The  seven  2.2  litre  jars  are  kept  in  place  by  a  steel  plate,  weighing 
approximately  32  kilogrammes. 

Each  sample  jar  has  a  water/sediment  intake  and  an  air  outlet. 

The  outlets  are  connected  to  a  common  outlet  which  ends  at  the 
top  of  the  glass  fibre  mast.  The  mast  also  helps  to  locate  the 
sandier.  Note  that  it  is  not  necessary  for  the  mein  outlet  to  be  out 
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of  the  water.  However, to  provide  sufficient  sampling  speed  it  should 
be  at  least  one  metre  above  the  top  of  the  jars,  and  the  main  out¬ 
let  must  be  kept  from  filling  with  water  before  the  sampling  begins. 

The  distance  from  the  intakes  to  the  nearest  part  of  the  sample 
jars  is  1.25m,  and  the  intakes  are  .65m  from  the  nearest  part  of  the 
sandier  that  touches  the  bed.  This  assures  a  minimum  of  disturbance 
of  the  bed  and  the  flow  near  the  intakes,  when  the  sampler  is  oriented  ' 

properly  relative  to  the  flow.  That  is,  the  bottles  are  placed  down¬ 
stream  from  the  intakes  if  a  net  current  is  present  and  the  horizontal 
symmetry  axis  of  the  sampler  is  aligned  along  the  wave  crests.  I 

The  most  common  working  depth  until  now  has  been  1.5m,  so  the 
driving  pressure  difference,  which  is  that  between  the  jar  stoppers 
and  the  end  of  the  air  outlet,  is  0.13  atmospheres.  Under  such  con-  1 

ditions  the  jars  fill  up  in  about  3.5  minutes,  which  corresponds  to 
an  intake  velocity  of  1.5m/s.  This  is  so  much  larger  than  the  typical 
settling  velocity  C'v-  .04o/s)  of  suspended  sand  that  one  can  expect  all  * 

grains  to  be  caught.  On  the  other  hand,  the  flow  velocities  induced 
by  the  suction  are  only  about  .03  m/s  at  positions  one  centimetre 
away  from  the  centre  of  the  intake,  so  the  suction  itself  causes  very  * 

little  disturbance  to  the  main  flow  field,  where  the  typical  velocity 
magnitude  is  0.5m/s. 

• 

TRANSDUCER  AND  PLOW  METER  CHARACTERISTICS. 

ftia  hydraulic  data  referred  to  in  thi3  paper  acre  obtained  from  • 

analysis  of  strip  charts , on  which  water  velocities  and  wave  heights 
were  recorded  during  the  sediment  sampling  process.  The  response 
time  of  the  chart  recorder  for  the  relevant  output  values  is  about  • 


80 


\ 

[ 

B 

I 

? 

g 

« 

I 

« 


0.2  seconds  and  thus  negligible  compared  to  the  time  scale  of 
natural  waves,  and  the  scaling  parameter  of  the  chart  recorder  is 
incorporated  in  the  calibration  of  the  instruments. 

Like  the  chart  recorder,  the  pressure  transducers  respond 
immediately  and  it  is  therefore  assumed  that  the  gain  is  the  3aue 
for  oscillatory  input  with  the  frequency  of  natural  waves  as  it  is 
for  static  pressures.  After  a  days  sampling  work r  the  transducer-* 
chart-recorder  system  is  calibrated  hydrostatically  before  dis¬ 
connection  in  order  to  achieve  maximum  accuracy  (±  .01m)  on 
determination  of  mean  water  level  and  wave  heights.  The  calibration 
arrangements,  is  shown,  in  Figure  4.2. 

The  flow  meters  have  been  calibrated  by  moving  them  through 
still  water  attached  to  a  heavy  pendulum  with  periods  between  2.5 
and  5.5  seconds.  This  calibration  3hows  that  the  instrument  res¬ 
ponses  are  linear  for  velocity  amplitudes  be twee en  0.1  and  2  m/s 
and  that  the  frequency  response  can  be  described  by 
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This  means  that  the  following  relation  exists  between  measured 

velocities  u  and  real  velocities  u: 
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Figure  4.2:  Instrument  pod  with  flowmeters  and  pre¬ 
ssure  transducer.  The  transducer  is  shown  with  the 
calibration  device.  Different  hydrostatic  pressure 
is  given  by  the  position  of  the  water  surface  in  the 
hose. 


The  equations  (4. 1)  and  (4.2)  provide  two  means  of  correcting  the 
velocity  measurements .  Equation  (4.1}  can  be  used  for  correcting 
spectra,  and  (4.2)  can  be  used  for  direct  correction  of  time  series 
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In  the  following  analysis ,  the  flow  meter  data  are  only 
used  to  determine  wave  period  and  time  average  velocities  U  and 
V  in  the  shore  normal  and  longshore  directions  respectively, 
while  the  oscillatory  velocity  amplitude  aw  is  obtained  from  the 
pressure  records  via  linear  wave  theory . 


1 


EXPERIMENTAL  PROCEDURE,  DETERMINATION  03?  c(z)  . 


a 


The  sampling  sites  are  chosen  so  that  the  hydraulic  conditions 
and  the  bed  topography  are  as  stationary  and  well  defined  as 
possible.  When  such  a  location  has  been  found  the  pod  with  flow 
meters  and  the  pressure  transducer  is  installed.  The  standard 
setup  has  two  ducted  impellor  flow  meters  and  one  pressure  trans¬ 
ducer.  The  shore  normal  flow  mater  is  placed  .2m  above  the  bed 
and  tiie  longshore  one  at  .3m.  The  pressure  transducer  is  normally 
placed  at  about  .5m  above  the  bed. 

The  sampler  is  placed  near  the  instrument  pod  so  that  the 
velocity  measurements  should  represent  the  flow  at  the  sampling 
spot  and  at  the  same  time  far  enough  away  (and  up  stream)  such  that  the 
sampler  and  the  pod  do  not  disturb  each  other. 

Then  the  position  of  the  intakes  is  carefully  adjusted,  and 
the  elevation  of  the  lowest  intake  above  the  bed  is  measured.  If 
the  bed  is  rippled,  the  intakes  are  always  placed  over  a  ripple 
crest. 

When  everything  is  ready,  a  clasp  which  has  been  blocking 
the  intake  and  outlet  hoses  is  removed,  the  sampling  starts.  When 
the  jars  are  full,  the  intake  elevation  is  measured  again  and  a  sand 
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sample  is  taken  from  the  bed  near  the  intakes.  If  the  bed  is  rippled, 
height  and  length  of  these  are  measured.  After  this  the  sampler  is 
brought  ashore  where  the  total  volume  of  the  samples  is  noted  and  the 
sand  is  filtered  out.  The  sand  is  dried  and  weighed  later  in  the 
laboratory  so  that  the  concentration  can  be  worked  out  from: 

-  _  mass  of  sand _  . .  ,, 

(density  of  sand)  x  (total  volume)  ' 


DETERMINATION  OF  THE  SEDIMENT  PARAMETERS:  d,  w  and  V 

The  bottom  sample  from  each  run  is  sieved  and  the  average  dia¬ 
meter  is  found  from 


Em.  d. 
r  i 

“5T 

i 


(4.4) 


where  m^  and  d^  are  mass  and  average  diameter  of  each  sieving  fraction. 
The  average  settling  velocity  is  similarly  found  from 


Em.  w. 
x  x 

Em, 


(4.5) 


Ae  variation  coefficient  v  of  the  settling  velocity  is  defined 


by 


V  -  Var{w}/D2 


(4.6) 


and  the  variance  is  found  from 


when  N  is  the  number  of  sieves.  The  settling  velocity  for  the 
different  sieve  fractions  is  determined  by  settling  tube  experi¬ 
ments  for  each  experiment  location.  This  is  necessary  because  the 
relation  between  grain  size  (d)  and  settling  velocity  (w)  varies  a 
great  deal  between  beaches,  see  Figure  4.3.  This  variability  is 
caused  by  differences  in  abundance  and  geometry  of  3hell  fragments 
in  the  bed  material.  These  shell  fragments  can  have  quite  large 
sieving  diameters  while  their  settling  velocity  is  much  smaller 
than  that  of  the  equivalent  quartz  sphere. 


0-2 


w  [m/s] 


Figure  4.3:  The  relation  between  sieving  grain  size  and  settling  velo¬ 
city  for  two  of  the  test  locations.  This  relation  varies  considerably 
from  beach  to  beach  due  to  variations  in  grain  geometry  and  mineral  con- 
position. 
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DETERMINATION  OF  h,  T,  H  ,  u  and  v. 

3 

Considering  the  scatter  of  the  concentration  measurements,  it 
was  concluded  that  a  representation  of  the  hydraulic  conditions  by  the 
water  depth  h,  the  significant  wave  height  the  average  period  T  and 
the  average  flow  velocities  u  and  V  would  be  detailed  enough.  The  wave 
period  is  defined  by  the  time  interval  between  successive  zero  down 
crossings  of  the  shore  normal  velocity  u (t) .  This  definition  is  chosen 
because  the  reversal  of  the  near  bed  velocity  is  essential  for  the  be¬ 
haviour  of  the  wave  boundary  layer. 

The  significant  wave  height  is  defined  as  the  average  of 
the  highest  third  of  the  zero  down  crossing  wave  heights.  That  is, the 
height  from  the  bottom  of  the  wave  trough  to  the  top  of  the  following 
crest.  The  water  depth,  h,  and  average  flow  velocities  u  and  v  are 
normally  found  by  running  the  pressure  transducer  and  flow  meters  in 
100  second  filter  mode  for  about  S  minutes  before  and  after  the  sampling. 
From  the  basic  parameters  mentioned  above  and  using  linear  wave  theory 
we  derive  the  semi  excursion  a: 

a 

a  .  - 3  , -  (4.3) 

2  sinh  —  h 

U 

and  the  velocity  amplitude  aw  using  u  ■  2ir/T.  The  mobility  number 
is  given  by 


I 


(aw) 2 
(s-l)  gd 


(4.9) 


where  s  is  the  relative  density  of  the  sand  and  g  is  the  acceleration 
due  to. gravity. 


0.05  a  *  n  c  0.25  a 


(4712) 


which  correaponds  to  a  diffusion  model  with  the  diffusivity  distribution 


'»<*>  -  S  * 


(4.14) 


via  the  one  dimensional ,  time  averaged  diffusion  equation 


dc 

e  — —  +•  WC 

a  dz 


(4. IS) 


£  is  an  apparent  diffusivity  because  equation  (4.14)  assumes  that  all 
suspended  sand  has  the  settling  velocity  w. 

The  best  fit  values  of  the  three  parameters  in  equation  (4.13) 
rewritten  as 
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(4.16) 
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are  found  by  the  least  squares  criterion 


EC  In.  c 


m  a:2 


0,  i  -  1,2,3 


(4.17) 
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2 

by  a  generalized  Newton  iteration  and  the  goodness  of  fit,  r  ,  is  de~ 
fined  by 


2  Cov  (in  3,  la  cl 
r  "  (Var  lln  c}  Var  l In  3))a's 


(4.18) 


If  the  bed  material  contains  fractions  with  different  settling  velocities 
and  if  all  grain  sizes  obey  the  diffusion  equation 
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(4.19) 


where  £  is  independent  of  grain  size,  the  concentration  of  different 
size  fractions  will  decrease  at  different  rates  away  from  the  bed 
and  the  shape  of  the  compound  concentration  profile  will  thus  de¬ 
pend  on  the  distribution  of  settling  velocities  for  the  bed 
material.  This  means  that  the  apparent  diffusivity  e^,  and 
in  particular,  which  comes  out  of  the  curve  fitting  process,  de¬ 
pends  on  the  w-distribution .  He  shall  now  show  how  this  effect 
can  be  eliminated  from  the  observations. 

In  the  following  we  assume  that  the  time  averaged  concentrations 

c(z,w)  of  sand  with  settling  velocity  w  obey  the  one  dimensional 

diffusion  equation  on  the  form  (4.19)  which  expresses  that  the 
dc  • 

upward  flux  equals  the  downward  flux  wc  at  any  level.  The 

general  solution  to  (4.19)  is 

z 

c(Z,w)  «  CQ(w)  exP  W  ^  ^*1^  (4.20) 

o 

where  CMw)  is  the  concentration  of  sand  with  settling  velocity  w 
at  z  •  0.  For  completely  homogeneous  sand  with  w  5  w  the  con¬ 
centration  profile  would  be  given  by 

z  - 

c(z)  -  C  exp  {-/ •  -jrdz,}  ,  w  =  w  (4.21) 

o 

For  convenience  we  introduce 
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/  fan 


(4.22) 


so  that  (4.21)  can  be  rewritten 


c(z)  »  C  e  for  w  =  w 
o 


(4.23) 


Natural  sediments  axe  not  homogeneous  and  the  variance  of 
settling  velocity  has  a  significant  effect  on  the  concentration 
profiles ,  making  the  gradients  larger  near  the  bed  and  smaller 
towards  the  surface. 

If  we  introduce  the  non  dimensional  settling  velocity 


w*  *  w/w 


(4.24) 


and  the  density  function  f(w*)  defined  by 


f(w')dw  ■  C  (w)/C  ,  0  <  w'  <  ® 

o  o 


(4.25) 


we  get  the  following  distribution  for  the  fraction  with  settling 
velocity  w 


c(z,w)  ■  C  a~W  <7f(w')dw' 

Q 


(4.26) 


and  for  the  total  concentration 


kz)  -  C  /  e“w'af  (w')dw’ 


(4.27) 


It  was  shown  by  Nielsen  (1979)  that  the  type  of  density  function 
is  not  critical  as  long  as  f(w')  has  the  right  variance  V  ■  Var(w)/w 
and  a  mean  value  of  1.  We  choose  a  r -distribution  be¬ 
cause  it.  fits  natural  w* -distributions  reasonably  well  (see 
Figure  4.5)  and  also  gives  simple  analytical  results 


f  (w* ) 


a/vr_ 

ra/v) 


exp  (-w*/V)w'  ^ 


under  the  assumption  that  £  and  thus  a  is  independent  of  w'. 
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The  desired  aim  is  to  find  a  method  for  determining  a 
true  E  from  a  measured  c-profile  of  non  homogeneous  sediment.  This 
can  be  done  in  the  following  way.  Equation  (4.22)  yields  the 
following  relation  between  ff  and  the  true  diffusivity 


w  do 

£  dz 


(4.30) 


do 

and  can  get  — 

QZ 


from  (4.29) : 


and  find 


(4.31) 


(4.32) 


Let  us  compare  this  true  diffusivity  with  the  apparent  one, 

e  ,  which  we  would  find  by  neglecting  the  variance  or  w.  The  apparent 
A- 

diffusivity  £  is  found  by  applying  (4.21)  directly  to  the  concen- 
trations  of  non  homogeneous  sediment  which  yields 


-  w 

d  In  8 
dz 


and  thus 


(4.33) 


(4.34) 
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The  result  shows  that  the  ratio  between  read,  and  apparent 
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diffusivities  at  a  given  level  depends  only  on  the  measurable  quantity 

c/Cq  and  on  the  variation  coefficient  Vy  which  is  very  convenient. 

We  see  that  the  effect  depends  strongly  on  V  as  one  would  expect  and 

that  the  effect  is  negligible  near  the  bed  (where  esc). 

o 

The  corrected  surface  diffusivity  Egc  is  derived  from  (4.14) 
and  (4.34) 
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DESCRIPTION  OF  EXPERIMENT  SITES 


PALM  BEACH 

Experiments  1  and  2  were  carried  out  at  the  southern  end  of 
Palm  Beach,  on  October  26,  1980.  Palm  Beach  is  the  most  northern 
of  Sydney's  beaches,  and  on  the  day  of  the  experiments  it  wets  in 
an  accreding  state  with  shallow  bars  moving  onshore.  The  bed  was 
covered  with  large  ripples  or  megaripples . 

CABLE  BEACH,  BROOME 

Experiments  3  through  24  were  performed  on  Cable  Beach  at  Broome, 
North  Western  Australia  in  November  1980.  Cable  Beads  has  a  tidal 
range  of  more  than  9  metres,  and  most  of  the  beach  profile  is  flat, 
without  bars.  Ripples  are  generally  small  or  absent  because  the 
sediment  is  very  fine.  Exceptions  from  this  general  pattern  are  some 
areas  near  the  neap  high  tide  line  where  the  sediment  is  coarse 
(d  *0.5  mm)  due  to  a  large  content  of  shell  fragments.  These  areas 
tend  to  develop  ridge  and  runnel  systems,  and  at  high  tide  the  bed 
is  covered  by  sharp  crested  vortex  ripples.  Experiments  9,  10  and  11 
were  performed  over  such  ripples.  The  morphology  and  wave  climate 
are  described  in  more  detail  in  Wright  et  al  (.1982) . 

MEROO  BEACH 

Experiments  25  through  35  were  carried  out  at  the  southern  end 
of  Meroo  Beach  approximately  8  kilometres  south  of  Ulladulla,  New 
South  Wales  in  February  1981.  The  beach  was  in  an  intermediate  state 
with  rather  confused  bar  topography.  The  bed  forms  were  large  vortex 
ripples  or  megaripples. 


SEVEN  MILE  BEACH 


Experiments  36  through  43  were  carried  out  at  the  northern  end 
of  Sevan  Mile  Beach,  south  of  Kiama,  New  South  Wales  in  March  1981 
and  experiments  61  through  66  were  done  at  the  same  location  in 
March  1982.  ®ie  northern  end  of  Seven  Mile  Beach  has  very  fine 
sand  throughout  and  is  exposed  to  the  prevailing  sea  and  swell  from 
the  Southeast.  The  modal  state  of  the  beach  is  therefore  the  extreme 
dissipative.  The  profile  is  flat,  bed-forms  are  normally  absent, 
and  the  sand  is  very  firmly  packed  (no  foot  prints) . 

EASTERN  BEACH 

Experiments  48  through  60  were  carried  out  on  Eastern  Beach  at 
Lakes  Entrance,  Victoria  in  May  1981.  At  the  time  of  the  experiments 
a  well  developed  bar  about  a  hundred  metres  from  the  shoreline  was 
moving  towards  the  shore.  The  depths  over  the  bar  crest  was  less  than 
half  a  metre  ah  low  tide  while  the  maximum  depth  in  the  trough  was 
about  two  metres.  All  suspension  measurements  were  taken  in  the  trough 
at  depths  of  approximately  one  metre.  The  bed  was  always  covered  by 
well  developed  vortex  ripples  that  showed  no  influence  at  all  from  the 
sometimes  strong  longshore  current  which  reversed  with  the  tide. 


WARRIWOOD  BEACH 

'  Experiments  67,  68  and  69  were  carried  out  in  April  1982  at 
Warriwood  Beach,  which  is  one  of  Sydney's  northern  beaches.  The  modal 
state  of  the  beach  is  intermediate  because  the  sand  is  quite  coarse 
and  the  headlands  at  both  ends  give  some  protection.  The  suspension 
measurements  were  taken  on  the  outer  slope  of  a  shallow  bar  and  the 
bed  was  covered  by  megaripples. 


WINDAWOPPA  BEACH 


Experiments  70  and  71  were  carried  out  on  Windawoppa  Beach  at 
Hawk's  Nest,  fifty  kilometres  northeast  of  Newcastle,  New  south  Wales 
in  July  1982.  'Hie  beach  was  recovering  after  being  severely  eroded 
by  storm  waves  the  preceding  week.  The  topography  was  rather 
confused  and  the  bed  was  covered  with  megaripples  at  the  test  site. 


PEARL  BEACH 

Experiments  44  through  47  were  carried  out  at  Pearl  Beach,  north 
of  Broken  Bay  near  Sydney  in  March  1981.  The  measurements  were  <-«ir«yn 
over  the  rippled  bed  about  10  metres  seaward  from  the  step  of  the 
reflective  beach.  Since  wave  reflection  is  very  significant  for 
locacion,  the  hydrodynamic  data  for  these  experiments  were  derived 
from  the  near  bed  water  velocities  rather  than  from  the  pressure 


variation. 
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The  vertical  length  scale  of  this  decay  (Zg)  is  approximately 
equal  to  the  ripple  height  n ,  so  the  concentrations  decay  by  a 
factor  100  ("  e4*5)  over  four  to  five  ripple  heights.  Figure  5.1 
shows  the  results  of  four  successive  runs  from  the  bar  trough  in 
Gippsland.  Apart  from  the  aboveman tioned  exponential  behaviour  the 
figure  also  shows  the  repeatability  of  the  tests. 

We  can  thus  describe  these  concentration  profiles  by  two 
parameters  Cq  and  £g.  Cq  determines  the  magnitude  and  determines 
the  distribution.  It  was  pointed  out  by  Nielsen  (1979)  and  it  may 
been  sen  in  the  data  of  Nakato  et.  al .in  Figure  5.2  that  concentration 
profiles  measured  over  the  ripple  trough  are  slightly  different  from 
those  over  the  crest.  Both  wilb.be  approximately  exponential  but  Cq 
will  be  1.5  times  smaller  over  the  trough  than  over  the  crest  while 
l  will  be  the  same  factor  larger r  so  that  the  total  amount  of 
suspended  material 


c 


cdz  • 


C  l 

o  s 


(5.1) 


is  the  same  for  all  vertical  sections.  All  the  field  measurements 
of  the  present  study  and  the  laboratory  experiments  quoted  in  the 
following  were  taken  over  the  ripple  crest. 

In  general  we  would  expect  l  to  be  a  function  of  the  bed  geometry 
the  boundary  layer  structure  and  the  parameters  that  determine  the 
sand  motion  in  the  boundary  layer  flow. 

When  the  bed  is  covered  with  sharp  crested  ripples,  the 
natural  vertical  length  scale  is  the  ripple  height  n. 


x  *  measurement  over  crest 

o  ■■  measurement  over  trough 
T  *  l'2  s  ,  a  »  0  051  m 
d  ■  0-14  mm 


c  ( m3/m3) 


Figure  5.2;  Tunnel  data  from  N aka to  at.  al.  (1977) . 

The  time-mean  concentration  decreases  faster  over  the  ripple 
crest  than  over  the  trough. 


The  boundary  layer  structure  will  depend  on  n /a  and  on  the  flow 
Reynolds  number  a2 u/v. 

The  flow  structure  is  dominated  by  the  lee  vortices  and  we 
found  in  Chapter  3  that  the  behaviour  of  a  sand  grain  in  vortex 
flow  is  determined  by  the  grain  Reynolds  number  wd/v  and  the 
velocity  ratio  RO/w.  The  observations  of  Tunsdall  and  Inman  (1975) 
show  that  the  velocity  scale  in  the  lee  vortices  is  au  so  we  may 
replace  Rfl  by  au. 


In  the  present  context  it  might  also  be  relevant  to  consider 
the  ratio  t^/T  between  the  spiralling  time  scale  (equation  3.51) 
and  the  wave  period. 

We  would  thus  expect  I.  to  be  given  by  a  relation  of  the  form 


The  following  empirical  considerations  shew  that  most  of  the 

variation  of  Z  /n  can  be  accounted  for  by  am/v  in  the  case  of  fully 
s 

developed  sharp  crested  ripples. 

Figure  S.3  shews  Z  /n  as  function  of  t  /T  and  we  see  that  for  all 

S  3 

the  data  where  at o/w  is  reasonably  big  (am/ w  >  9)  the  variation  i3  not 
significant  except  for  a  few  data  points  with  t^/T  >  1.  The  points 
that  fall  far  below  correspond  to  small  values  of  am/ w  (  <  7.5)  which 
indicates  that  for  am/m  less  than  about  9  the  vortex  trapping 
mechanism  is  less  efficient.  TO  calculate  tg  we  have  used  equation 
(3.51)  and  assumed  that  the  velocity  scale  in  the  vortex  is  am  and 
that  the  vortex  radius  is  0.5n.  We  thus  find  the  angular  velocity  8 
in  the  vortex  to  be 


8 


a  2ir 

0.5n  T 


and 
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(5.3) 


(5.4) 


Experiments  with  extremely  large  values  of  t  /T  will  also 

s 

lead  to  unusually  high  values  of  Z  /r\.  Large  values  of  t  /T  occur 

3  S 

when  t  and  w  are  both  small  and  n /a  is  large.  In  relative  terms, 
such  conditions  are  characterized  by  small  grain  Reynolds  numbers 
(wd/v  s  1)  and  small  flow  Reynolds  numbers  (a2m/v  $  2000) , 


e.g.  the  experiments  "H"  and  "N"  in  Figures  5.3  through  5.5. 
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Figure  5.3s  The  variation  of  l3/rt  with  ts/T  is  insignificant  for 
t  /T  <1,  which  includes  all  field  observations.  But  the  variation 
wfth  au/w  must  be  strong  since  all  data  with  aw/w  <  7.5  fall 
significantly  below  the  rest.  Since  the  typical  value  of  l%/ n  i* 

1.4  a  l/ln2 ,  we  can  derive  the  following  rule  of  thumb: 

The  concentration  decreases  by  a  factor  two  when  the 
elevation  is  increased  by  one  ripple  height. 

The  symbols  used  in  Figure  5.3  and  the  following  are  defined  in  Table 

5.1.  The  field  data  are  primarily  described  by  the  bed  topography: 

R  for  vortex  ripples,  T  for  small  ripples  with  low  steepness  due  to 

high  shear  stresses,  F  for  flat  beds  and  MR  for  megaripples.  But 

they  may  alternatively  be  described  by  the  wave  characteristics:  N 

for  non  breaking,  S  for  spilling  breakers  or  P  for  plunging  breakers. 


Synbol  Author  Have  Period  Grain  Size  Facility 

(Seconds)  (nun) 
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Figure  5.4;  The  relative  vertical  suspension  length  scale,  J3/n,  as 

function  of  the  velocity  ratio  aw/ w.  I  /n  grows  with  aw/ w  between  3 

and  10  as  the  vortex  trapping  becomes  mors  efficient.  I  /n  becomes 

constant  for  au/w  >  10,  probably  because  the  ripples  get3more  rounded. 

If  the  ripple  "sharpness"  is  maintained  as  in  MacDonalds  experiments 

(•),  l  /r\  may  continue  to  grow, 
s 

When  the  ripple  height  is  snail  (n/<z  S  0.01)  as  for  the  data  marked 
(T)  l  /n  will  be  luge  because  l  does  not  become  substantially  smaller 
than  the  total  bed  roughness.  9 

Megaripples  (MS)  do  not  shed  strong  vortices,  so  the  concentration 
profiles  do  not  scale  on  the  megaripple  height  but  on  the  flat  bed 
roughness . 
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Figure  5.4  shows  £g/n  versus  aw/w.  The  main  trend  here  is  that 

Zg/ri  increases  with  am/ w  for  3  <  am/ w  <  10  ,  indicating  that  the  vortex 

trapping  mechanism  gets  more  effective  with  increasing  velocity 

ratio.  For  am/ w  >  20r £  /n  tends  to  decrease  again,  probably  because 

the  ripples  get  more  rounded  and  less  efficient  vortex  makers.  The 

dots  towards  the  upper  right  hand  corner  are  measurements  from 

MacDonald  (1977)  who  used  solid  half  circular  "ripples"  with  just  a 

small  amount  of  loose  sand  around.  These  wooden  "ripples"  will  of 

course  maintain  their  shape  and  vortex  making  capability  independent 

of  am/ w  so  that  l  /n  can  continue  its  upward  trend.  (The  shown 
s 

"MacDonald  data”  all  have  a/n  >  14.  When  a  is  smaller  (4.9  <  a/ n  <  13) 

it  seems  to  restrict  the  increase  of  l  /n  somewhat) . 

3 

The  data  marked  "T"  represents  conditions  where  the  relative 
ripple  height  (n /a)  has  decreased  due  to  increasing  <|i  and  9* ,  and 
become  smaller  than  the  equivalent  flat-bed-roughness,  r.  Undt ■ 
such  conditions  igA i  is  again  large  because  Zg  now  scales  on  r. 

The  influence  of  t  /T  on  I  /r\  is  shown  by  the  circled  N-experiment 
in  Figure  5.4  but  possibly  also  by  the  experiments  of  Nakato  et.  al. 
which  are  marked  I.  Those  experiments  were  performed  with  constant 
au/w  and  three  different  periods:  1.2,  1.3  and  2.4  seconds.  They 
show  a  considerable  increase  of  l  /n  with  T  which  may  be  due  to  the 

3 

related  increase  in  t  /T  -  gT/w. 
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Figure  5.5:  Observed  relative  ripple  heights  from  the  present 
field  study  and  from  some  of  the  laboratory  studies  quoted.  The 
curve  represents  the  trend  of  all  available  laboratory  experiments 
with  quartz  sand  (Nielsen  1981) .  The  bed  forms  referred  to  as  I 

megaripples  have  rounded  crests  and  are  not  shedding  vortices 
periodically.  Extreme  values  of  n/a  (>0.25)  tend  to  occur  for 
small  grain  -  and  flow  Reynolds  numbers . 
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GRAIN  SIZE  VARIATION  WITH  ELEVATION  OVER  VORTEX  RIPPLES 


The  profiles  measured  in  tests  57  through  60  axe  very  similar 

(see  Figure  5.1),  and  it  was  thought  reasonable  to  pool  saiqples  from 

the  same  elevation  to  get  amounts  big  enough  for  a  sieving  analysis. 

The  result  of  the  sieving  is  shown  in  Table  5.2.  The  distributions 

from  the  four  highest  elevations  are  so  similar  that  they  could  not 

be  distinguished  in  a  usual  cumulative-distribution-plot.  The 

distributions  from  the  lower  levels  differ  somewhat  more,  mainly 

with  respect  to  the  coarsest  material  (d  >  0.5  mm) ,  but  the  general 

picture  is  that  the  typical  grain  size  (  and  settling  velocity)  varies 

much  less  with  elevation  than  a  diffusion  modal  would  suggest. 

If  the  entrainment  process  was  diffusion  with  s(z)  independent 

of  grain  size  and  if  settling  velocities  were  r-distributed  with 

variation  coefficient  V  at  r  .  Then,  using  the  terminology 

o 

introduced  in  'action  4,  we  find  that  the  mean  settling  velocity 
will  vary  with  z  as 


w(z) 


w(z  ) 
o 

1-KTV 


(5.5) 


where  a  is  the  non  dimensional  elevation 


z 

f  W( 
J  e 


V 


e(z) 


dz 


(4.22) 


The  variation  corresponding  to  (5.5)  with  z  ■  0.01  m,  w(z  )  *  0.053m/s, 

o  o 

V  »  0.106  and  e  *  0.0043  tor/s  is  shown  in  Figure  5.6  together  with  the 
observed  values  from  Table  5.2.  We  see  that  the  observed  values  aura 
much  more  constant  than  equation  (5.5)  suggests,  which  indicates  that 


Table  5.2 t  Results  of  sieving  analysis 
bed  and  from  the  seven  different  intakes 
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Figure  5.6:  Predicted  decrease  of  w(z)  from  diffusion 
model  compared  to  the  observed  values  (x) .  The  observa¬ 
tions  show  that  9(z)  is  much  more  constant  than  the  dif¬ 
fusion  model  predicts. 


no 


the  distribution  of  individual  size  fractions  are  less  dependent 
of  settling  velocity  than  the  diffusion  model  predicts 
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(5.6) 


Another  way  of  illustrating  this  is  by  plotting  concentration  profiles 
for  individual  size  fractions  and  see  whether  they  have  different 
slopes  in  accordance  with  (5.6).  Ibis  has  been  done  in  Figure  5.7 
and  we  see  that  the  differences  in  slope,  especially  above  z  -  0.08m 
are  much  smaller  than  predicted  by  (5.6).  Zn  fact  all  the  slopes  are 
identical  for  z  >  0.8,  indicating  that  all  grain  sizes  are  distributed 
alike,  irrespective  of  settling  velocity.  The  mechanism  responsible 
for  this  cannot  be  diffusion  but  must  be  a  convective  process  where 
the  sand  grains  travel  with  the  released  lee  vortices  after  being 
trapped  as  explained  in  Chapter  3.  The  effect  described  above  has 
been  noticed  before  for  exasple  by  Coleman  (1970)  who  studied 
concentration  profiles  for  different  sand  sizes  in  steady  flow. 
Maintaining  the  diffusion  terminology  he  concluded  that  big  grains 
experience  a  much  larger  diffusivity  than  small  ones  in  the  same  flow. 
It  is  more  profitable  however  to  discard  the  diffusion  terminology  when 
the  flow  has  a  pronounced  vortex  structure,  like  over  ripples  and 
behind  dunes,  and  describe  the  concentration  profile  in  terms  of  the 
length  scale 

v*>  -  E  s_aK[  (5-7) 


which  is  very  similar  to  the  corresponding  length  scale  for  vertical 
velocity  fluctuations 
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as  indicated  by  the  measurements  of  MacDonald  (1977)  and  Nakato 
et.  al.  (1977),  see  Figures  2.6  through  2.8. 

Grain  size  variation  with  the  distance  from  the  bed  under 
breaking  waves  w as  measured  by  Kana  (1978)  who  found  a  more  pronounced 
variation  in  average  grain  size  between  10  and  60  centimetres  above 
flat  beds.  This  shows  that  the  entrainment  mechanism  is  different 
outside  the  wave  boundary  layer,  under  breaking  waves,  and  probably 
more  like  diffusion. 


CONCENTRATION  PROFILES  OVER  FLAT  BEDS 

Flat  beds  of  firmly  packed  sand  axe  common  in  the  surf  zone 
of  flat,  dissipative  beaches. 

under  such  conditions  the  concentration  gradients  are  very  large 
near  the  bed.  Typicaly  the  concentrations  will  decrease  by  am  order 
of  magnitude  over  the  lowest  ten  centimetres. 

Figure  5.8  shows  two  profiles  measured  under  very  similar 
conditions  at  Seven  Mile  Beach,  March  1982. 

The  suspended  sand  near  a  flat  bed  forms  long  stream  parallel 
clouds  twice  every  period  near  the  velocity  extremum  in  either 
direction;  more  dense  under  the  stronger  shoreward  flow.  These 
streaky  clouds  break  up  when  the  flow  reverses  and  some  of  the  sand 
travels  to  slightly  higher  elevations  with  the  turbulent  water,  but 
most  of  the  sand  will  have  settled  out  before  the  next  velocity 
extremum  occurs  and  forms  new  clouds. 
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CONCENTRATION 

Fiqure  5.8:  Concentations  measured  over  a  flat  bad  under  plunging 
BSSBErar- The  curve  is  the  result  of  the  fitting  procedure  described 

in  connection  with  equation  4.13.  C  »  1.08  x  10“2,  £  “  1.86  x  10"3  m2/s 

£s  «  4.0  x  10“2  m2/s,  r2  *  0.954. 
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The  width  and  thickness  of  the  clouds  is  normally  3-5  centimetres 
under  field  conditions  and  results  in  vertical  length  scales  of  the 
same  magnitude  near  the  bed. 

The  structure  of  the  sand  clouds  looks  very  similar  to  that  of  the 
"sublayer  streaks"  in  steady  flow.  See  Cantwell  (1981) .  The  characteristic 
£g  value  of  3-5  centimetres  is  approximately  equal  to  the  hydraulic 
roughness,  r  given  by  equation  (2.65)  or  about  a  hundred  grain  diametres. 

The  length  scale  will  normally  increase  quite  rapidly  with  distance 
from  the  bed.  under  breaking  waves,  the  increase  is  due  to  mixing 
induced  by  wave  breaking  and  concentrations  of  the  order  of  magnitude 
10“ 4  may  extend  all  the  way  to  the  surface.  Under  non  breaking  waves 
the  steepening  of  the  profiles  occurs  at  much  lower  concentations  (of 
the  order  of  10“5)  which  are  made  up  mainly  by  the  very  finest  fractions 
of  the  bed  material,  see  Figure  5.9. 
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Figure  5.9:  Sediment  concentrations  over  a  flat  bed  under  non 
breaking  waves,  under  non  breaking  waves,  the  concentration  magnitude 
will  drop  to  somewhere  between  10“  *♦  and  10“ 5  before  the  profile  gets 
steeper. 


Figure  5.10  shows  7-  /d  (»  s  /wd)  for  all  measurements  carried 
out  over  flat  beds  or  megaripples. 

He  see  that  the  near  bed  length  scale  determined  from  the 
curve  fitting  as  e  /»,  is  generally  about  a  hundred  grain  diameters 

B 

or  of  the  same  magnitude  as  the  hydraulic  roughness,  r,  of  a  loose 
sand  bed  as  derived  from  the  friction  measurements  of  Cars Ians  et.  al. 
(1969)  and  Lofquist  (1980) 


r  -  8  n 2/X  +  190  /e*  -  0.05'  d 


(2.65) 


which  gives 


r/d  -  190  /9'  -  0.05' 


for  a  flat  bed. 
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Eigure  5.10:  The  value  of  l3  near  the  flat  bed  is  about  a  hundred 
grain  diametres  or  approximately  equal  to  the  hydraulic  roughness 
given  by  equation  (2.65).  There  is  a  slight  tendency  for  plunging- 
breakers  (P)  to  give  larger  near  bed  length  scales  than  spilling 
breakers  (S)  indicating  that  the  turbulence  of  a  plunging  breaker 
penetrates  further  towards  the  bed. 

The  data  in  Figure  5.10  show  no  systematic  trend  and  the 

; 

>  relation  between  l%  and  r  may  be  expressed 


Ijx  -  0.81  t  0.50 


(5. 1C 
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the  standard  deviation  is  very  large  indeed  and  this  is  of 
course  partly  due  to  the  fact  that  the  s angler  is  not  really  designed 
for  sampling  over  flat  beds.  There  are  too  few  intakes  within  the 
lowest  ten  centimetres,  and  it  is  very  difficult  to  adjust  and 
measure  the  elevation  of  the  lowest  intake  with  sufficient  accuracy. 
There  are  however  no  other  detailed  measurements  available  at  the 
moment  as  far  as  the  author  knows. 

One  might  have  expected  the  appropriate  length  scale  to  be  the 
boundary  layer  thickness  or  the  displacement  thickness  5  ^  rather  than 
the  rou^mess  r,  and  indeed  Zg/<5 1  is  typically  of  the  order  of 
magnitude  one. 

However  Zg  does  not  seem  to  grow  with  a  in  the  way  S i  does. 

(For  the  relevant  range  of  r/a  we  have  5 1  a'75  r‘2S) .  Figure  S.ll 

shews  Ig/61  versus  a/d  for  the  same  data  as  plotted  in  Figure  5.10. 

The  scatter  is  very  large  indeed,  but  we  see  the  trend  of 
Zg/$1  to  decrease  with  increading  a/d,  which  probably  means  that 
l  is  fairly  independent  of  a  while  grows  like  a  to  a  power 
between  0.25  and  0.75.  (For  the  determination  of  Sj,  see 
equations  3.61,  2.62  and  2.65) . 
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Figure  5.11;  Relative  near  bad  length  scales,  Zs/5i  as  function  of 
a/ d  from  measurements  over  flat  beds  or  megaripples.  The  height 
n  of  the  megaripples  Is  not  Included  In  the  calculation  of  S\. 


CONCENTRATION  PROFILES  OVER  MEGARIPPLES 

Wh en  the  bed  forms  axe  nicely  regular  like  long  crested  vortex 
ripples  the  suspended  sediment  pattern  is  correspondingly  regular, 
uniform  in  the  longshore  direction  and  periodic  in  the  shore-normal 
direction.  It  is  however  quite  common  to  find  different,  more 
irregular  bed  forms  like  megaripples  and  the  suspended  sediment 
pattern  over  them  shows  strong  and  complicated  spatial  and  temporal 
variability. 

Megaripples  are  irregular  bedforms  with  typical  heights  of  ten 
to  thirty  centimetres  and  lengths  of  one  to  two  metres  in  the  shore 
normal  direction.  The  crest  lengths  are  of  the  same  order  of  magni¬ 
tude  as  the  crest  to  crest  length,  so  the  pattern  is  three  dimensional 
and  quite  irregular. 

The  crests  are  very  rounded  compared  to  those  of  fully  developed 
vortex  ripples  and  there  is  no  periodic  formation  of  lee  vortices. 

The  suspension  pattern  is  therefore  in  most  pleaces  rather  like  that 
over  a  flat  bed  with  sand  entrained  by  turbulence  bursts  that  form 
long  streaks  parallel  to  the  wave  motion  and  reaching  heights 
comparable  to  the  roughness  length,  r,  of  a  flat  sand  bed.  However 
in  some  places  and  under  some  waves  large  sand  fountains  are  formed 
by  strong  vortices  that  are  formed  around  the  current  "myi™  and 
released  by  the  following  flow  reversal  to  travel  upward  carrying 
suspended  sand  up  to  elevations  as  high  as  one  or  two  metres  above 
the  bed.  Horizontal  sections  through  the  fountains  will  typically 
have  diameters  of  0.3  to  0.5  metres. 

Megaripples  are  often  non  stationary  with  a  continually 
changing  pattern  so  the  distribution  of  active  (in  the  sense  of 
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Figure  5.13  shows  time  series  of  shore  normal  velocity  at 
z  =»  0.20m  together  with  instantaneous  sediment  concentrations 
measured  at  z  »  0.04m  near  the  sampling  location  of  the  two  c(z) 
profiles  shown  in  Figure  5.12.  During  the  first  three  minutes  of 
the  record,  the  detector  was  at  a  very  active  spot  where  strong 
bursts  of  suspended  sand  were  created  by  every  second  or  third 
wave  and  then  decayed  over  a  twenty  to  thirty  second  period. 

Over  the  last  3  minutes  nothing  much  happened  in  the  way  of 
sdeiment  suspension  although  the  waves  were  unchanged.  The  active 
spot  had  died  or  moved  on. 


,^u(a2o,tl(m/sI 


rilWIIJ 

HU 


mam 


<(004,t)  BY  VOLUME 


h'’  :  i ' i i i  ] 

Figure  5.13:  Simultaneous  time  series  of  shore  normal  velocities  and 
instantaneous  sediment  concentrations.  For  the  first  3  minutes,  the 
area  around  the  sediment  detector  was  much  more  active  than  during 
the  last  half  of  the  record.  There  is  no  corresponding  change  in  the 
water  motion,  so  the  change  must  be  due  to  changes  in  the  megaripple 
topography.  The  instantaneous  concentrations  were  measured  via  gamma 
radiation  absorption  by  Alan  Davison  of  the  Australia  Atomic  Energy 
Commission . 
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Hie  near  bed  sediment  distribution  over  megaripples  away  from 
the  active  spots  is  similar  to  that  over  a  flat  bed.  The  shape  of  the 
upper  part  of  the  profiles  is  similar  to  that  under  breaking  waves 
over  flat  beds  and  the  present  data  does  not  enable  us  to  distinguish 
the  effects  of  breaking  waves  from  those  of  sand  fountains . 


THE  DIFFUSION  APPROACH  TO  NEAR  BED  SEDIMENT  DISTRIBUTIONS 

The  classical  approach  to  suspended  sediment  phenomena  is  the 
gradient  diffusion  approach  where  the  upward  sediment  flux  is 
assumed  proportional  to  the  concentration  gradient  ~  and  to  the 

3Z 

diffusivity  e,  which  may  be  a  function  of  z  but  should  be  essentially 
independent  of  w  and  approximately  equal  to  the  eddy  viscosity, 

T 

if  diffusion  of  sediment  is  physically  analogous  to  diffusion  of 
momentum. 

The  distribution ' s  dependence  on  w  enters  through  the  assumption 
of  local  equilibrium  between  the  upward  flux  e  t— •  and  the  downward 

OZ 

flux  wc.  This  equilibrium  is  expressed  by  the  time  averaged 
diffusion  equation 


+  wc 


(4.19) 


It  should  be  noted  here  that  Nakato  et.  al  (1977)  found  that 
there  is  no  local  equilibrium  between  time  averaged  upward  and 
downward  fluxes  over  ripples.  The  net  flux  is  upward  over  ripple 


crests  and  downward  over  the  troughs.  However  the  diffusion  equation 
may  still  be  applied  in  a  space  averaged  sense.  That  is  if  c  is 
defined  as  the  average  over  a  ripple  length. 

If  the  diffusivity  is  essentially  the  same  as  the  eddy  viscosity 
it  should  be  scaled  by  a  typical  velcoity  vv  and  a  typical  length 
lv  of  the  boundary  layer: 


v  v 


const 


(5.11) 


Nielsen  (1979) tested  several  relations  of  this  type  for  the 
case  of  rippled  beds  with  negative  conclusions.  The  diffusion  model 
does  not  apply  over  rippled  beds  because  the  dominating  entrainment 
mechanism  is  convection  through  the  entire  boundary  layer  thickness 
of  sand  trapped  in  vortices.  We  have  studied  this  process  in 
previous  sections. 

The  turbulence  structure  over  flat  beds  is  evidently  different, 
so  it  is  not  clear  a  priori  whether  the  entrainment  process  is  mainly 
diffusion  or  convection.  That  is  whether  the  distributions  will  be 
determined  by  a  diffusivity,  given  by  the  turbulence  structure 

e  «  v  (5.12) 

T 

leading  to  different  length  scales  e/w  for  different  grain  sizes, 
or  the  distributions  of  all  grain  sizes  will  have  a  common  vertical 
length  scale,  l  ,  determined  by  the  boundary  layer  structure. 
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We  have  found  in  the  previous  section  that  a  common  length  scale 
namely  the  bed  roughness  applies  pretty  well: 


l 


0.81  ±  0.50 


(5.10) 


in  accordance  with  a  convection  model. 

However  the  scatter  is  very  considerable  and  in  fact  a 
diffusion  model,  assuming 


constant 


(5.13) 


I 


does  not  fall  through  completely  on  the  basis  of  the  present  data.  I 

Figure  5.14  shows  sa/u#5 i  versus  w  for  the  same  megaripple  -  (MR) 
and  flat  bed  data  (F)  that  were  used  in  Figure  5.10.  We  see  that  the 
typical  magnitude  of  eQ/u#<Sj  is  around  0 . 3  which  we  would  expect  in  | 

view  of  Jons son  and  Carlsen's  eddy  viscosity  measurements  that  aura 
shown  in  Figure  2.10.  However,  there  is  a  significant  tendency  for 
the  observed  eddy  viscosities  to  increse  with  w.  This  is  evidence  I 

in  favour  of  the  convection  process  because  the  universal  (<■  independent 
of  w)  length  scale  l  of  the  convection  model  corresponds  to  e 

S  B 

being  proportional  to  w: 


w  l 


(5.14) 
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Figure  5.14:  Measured  relative  diffusivlties  for  flat  bed  (F) 
and  mega ripple  (MR)  data.  The  typical  magnitude  is  0.3  which 
one  would  expect  from  Jonsson  and  Carlsen's  eddy  viscosity 
measurements ;  but  the  observed  diffusivities  increase  with 
w,  which  indicates  that  the  entrainment  process  is  reaily 
convective . 


CONCENTRATION  MAGNITUDE 


The  concentration  magnitude  is  determined  by  CQ  which  is 
defined  as  the  time  averaged  concentration  (by  volume)  at  the 
ripple  crest  level  and  it  is  determined  experimentally  by 
extrapolation  of  the  fitted  profiles  given  by  4.10  or  4.13  to 
the  level  z  ■  0. 

It  is  natural  to  expect  CQ  to  depend  mainly  on  the  ratio 
between  moving  -  («  shear  stress)  and  stabilizing  forces  on  the 
bed  sediment  particles. 

Thi3  ratio  is  expressed  by  the  Shields  parameter 
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the  form  drag  on  ripples  and  other  bedforms  is  unlikely  to  influence 
the  motion  of  individual  grains  so  roughness  due  to  bed  form  geometry 
i.e.  the  first  term  of  equation  2.65,  should  not  be  included  in  the 
calculation  of  the  Shields  parameter.  The  roughness  due  to  moving 
grains  is  mainly  a  function  of  8'  as  we  found  in  connection  with  (2.65) 
so  any  effect  of  moving-grain-roughness  is  automatically  included  if 
we  assume  a  relation  of  the  following  form 


C  -  F(8*) 
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where 
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and  S'  is  calculated  as  f  from  (2.62)  with  r  •  2.5d 
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Figure  5 . 15  shows  Cq  as  function  of  9 '  for  all  available 
laboratory  data.  The  concentrations  measured  by  Homma  et.  al. ,  (J) , 
are  generally  three  to  four  times  larger  than  the  main  line  of  the  rest 
which  remains  unexplained.  They  were  measured  by  the  same  kind  of 
optical  device  as  used  by  Nakato  et.  al.  (I) .  The  very  different 
grain  diameters  (0.082  -  0.55  mm)  used  by  Nielsen  (1979)  follow  the 
same  trend,  so  the  variation  of  relative  grain  size  (a/ d)  is 
properly  accounted  for.  However  there  is  a  significant  tendency 
for  data  with  longer  periods  or  higher  flow  Reynolds  numbers  (a2u/v) 
to  fall  below  those  with  smaller  Reynolds  numbers .  See  for  example 
the  data  of  Nakato  et.  al.  that  are  marked  "I".  The  three  experi¬ 
ments  were  carried  out  with  the  same  sand  and  the  same  velocity 
amplitude  (am)  but  with  different  periods:  1.2,  1.8,  and  2.4  seconds 
and  correspondingly  different  Reynolds  numbers:  1-4*104,  2.1*104,  and 
2.7*  104 . 

The  power  function 

C  «  0.036  0a»7  (5.18) 

o 

is  found  by  log- log  linear  regression  and  redrawn  in  Figure  5.16 
to  illustrate  the  fact  that  all  the  field  data  fall  significantly 
below  the  main  stream  of  the  laboratory  data.  The  field  data 
correspond  to  much  longer  periods  than  the  laboratory  data,  typically 
8  seconds  as  opposed  to  2  seconds,  and  will  therefore  have  much  larger 
Reynolds  numbers  for  the  same  9 ' .  The  discrepancy  between  experiments 
with  different  Reynolds  nuafcers  is  caused  by  the  way  in  which  9*  is 
calculated.  The  use  of  equation  (2.62)  assumes  rough  turbulent  flow 


Figure  5.16:  C  versus  9'  for  fiald  data.  The  trend  is  the  same  for 
ripples,  aegari§ples  and  flat  beds  (R,  MR  and  F)  but  all  field  data 
fall  significantly  below  the  line  given  by  (5  ..18) ,  which  represents 
the  general  trend  of  the  laboratory  data. 


conditions,  while  viscous  effects  are  probably  still  quite  important 
in  most  of  the  laboratory  experiments .  The  viscous  contribution  to  the 
friction  factor  may  be  roughly  estimated  by 
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w 


(2.53) 


which  we  derived  above  for  laminar  flow  over  a  smooth  bed.  This  will 
be  considerably  larger  than  the  values  derived  from  (2.62)  for  small 
scale  experiments  and  therefore  the  proper  Shields  parameter  for  the 
laboratory  experiments  is  larger  than  9 • .  In  Figure  5 . 17  the  laboratory 
data,  have  been  moved  to  the  right,-  to  closer  agreement  with  the  field 
data  by  including  viscous  effects  in  the  friction  factor  and  the 
Shields  parameter  in  the  following  way.  For  the  friction  factor  we 
use  the  combined  formula 
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(5.19) 


and  get  the  revised  Shields  parameter 


9*  “  0.5  ff#*  (5.20) 

The  Reynolds  numbers  for  field  data  range  roughly  between  10s  and 

106  so  in  their  case  9#  is  practically  identical  to  9'.  The  scatter 

is  still  very  large,  especially  towards  the  lower  values  of  C  and 

o 

unfortunately  it  is  real.  As  mentioned  in  connection  with  Figure  5.12 
the  concentrations  can  be  an  order  of  magnitude  different  due  to 
details  in  the  bed  form  topography  under  the  same  waves  and  with 
the  same  bed  sediment,  so  if  we  are  looking  at  general  relations 


FI  crura  5.17;  C  as  function  of  8t  for  field  and  laboratory  data. 
The  major  deviations  from  the  main  trend  axe  found  for  the 
laboratory  measurements  of  Homma  at.  al  (J)  and  for  field 
measurements  over  hard,  flat  bads  at  moderate  values  of  the 
Shields  parameter  (40,  41,  42  and  65) . 
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between  Sheilds  parameter  and  concentration  magnitude  we  must  expect 
this  kind  of  scatter  for  beds  covered  with  irregular  bed  forms. 

For  flat  beds  and  rippled  beds  with  samples  consistently  taken 

over  the  ripple  crest,  we  might  expect  a  somewhat  smaller  scatter 

and  that  is  generally  what  we  find.  However  over  soma  flat  beds  of 

fine  sand  we  tend  to  find  concentrations  that  are  extremely  small 

compared  to  the  general  trend.  Such  concentrations  were  found  in 

experiments  40,  41,  42  and  65  at  Seven  Mile  Beach.  The  fine  sand 

(d  3  0.17  mm)  was  very  firm  to  stand  on,  nearly  like  solid  concrete, 

the  bed  was  completely  flat  and  nearly  no  sediment  motion  occurred 

although  the  Shields  parameter  was  larger  than  0.55.  That  is  an  order 

of  magnitude  larger  than  the  conventional  critical  value  for  the  onset 

of  sediment  motion  (8  3  0.05) . 

c 


INFLUENCE  OF  WAVE  BREAKING  ON  CONCENTRATION  PROFILES 


The  wave  boundary  layer  turbulence  decays  very  rapidly 
with  distance  from  the  bed  as  described  by  Tunsdall  and  Inman 
(1975)  and  measured  by  MacDonald  (1977)  and  Nakato  et.  al.  (1977) , 
see  Figure  2.19.  The  data  support  the  following  rough  estimate 
of  the  scale  of  turbulent  velocities 


0.9  u# 
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(2.81) 


This  means  that  the  suspended  sediment  distribution  outside  the 
boundary  layer  (z  2  56 i)  is  likely  to  be  determined  by  entrainment 
and  mixing  due  to  sources  such  as  wave  breaking. 

Different  types  of  breaking  waves  are  likely  to  generate  mixing 
of  vastly  different  character  and  strength.  Waves  that  plunge  heavily 
on  shallow  barm  or  on  the  step  of  steep  beaches  can  form  very  strong 
jets,  that  penetrate  right  through  to  the  bed  and  thus  introduce 
very  strong  external  turbulence  into  the  boundary  layer  itself. 

The  jet  is  also  able  to  inject  large  amounts  of  entrained  air 
into  the  boundary  layer,  and  when  this  air  rises,  it  generates  large 
localized,  upward  water  velocities  that  act  as  very  efficient  elevators 
for  suspended  sand.  The  magnitude  of  the  vertical  water  velocities, 
generated  in  this  way  can  be  estimated  from  the  typical  heights,  Z, 
reached  by  the  water  and  sand  dragged  up  by  the  escaping  air. 

It  is  not  uncommon  to  see  splashes  of  white  or  brownish  water  rise 
about  a  metre  above  the  local  water  ievel  shortly  after  a  wave  has 


plunged 
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The  relation  between  starting  velocity  v^  and  culmination 


height  Z  is 


(5.21) 


so  even  with  a  conservative  estimate  of  Z  »  0.5  m  we  find 


v  3  3  m/s 

o 


(5.22) 


which  is  a  hundred  times  the  typical  sediment  settling  velocities. 

These  large  velocities  of  course  are  concentrated  around  large  escaping 
air  pockets  and  of  short  duration  but  they  dominate  the  sand  entrain¬ 
ment  in  the  areas  around  the  plunge  point.  The  described  process 
may  create  some  very  strange  looking  concentration  profiles  when  the 
rising  plumes  spread  along  the  surface  above  areas  with  little 
sediment  motion.  Then  one  may  observe  that  the  concentrations 
increase  toward  the  star  face. 

Mixing  due  to  spilling  breakers  and  bores  has  a  very  different 
character.  The  writer  once  observed  the  effect  of  spilling  breakers  in 
a  wave  tank  (T  •  1.7s,  h  »  0.40  m,  H  *  0.2  m,  d  3.082  urn)  at  the 
Institute  of  Hydrodynamics  and  Hydraulic  Engineering  (ISVA) ,  Technical 
University  of  Denmark,  during  the  following  experiment:  At  first 
the  waves  were  not  breaking  and  the  suspended  sand  formed  a  well  defined 
layer  over  the  rippled  bed,  only  a  few  centimetres  thick.  Then  the 
aavlituda  of  the  wave  generator  was  slightly  increased  so  that  the 
waves  began  to  break  and  passed  over  the  test  section  as  gently 
spilling  breakers,  the  wave  height  being  practically  unchanged. 
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The  suspension  distribution  was  unchanged  for  four  or  * 

five  wave  periods  after  the  onset  of  wave  breaking.  Then  it  changed 
very  rapidly  and  apparently  reached  a  new  equilibrium  with  con¬ 
siderable  concentrations  throughout  the  water  column.  The  transition  t 

took  only  about  two  wave  periods  (from  a  visual  judgement). 

When  the  wave  maker  was  turned  down  again  so  that  the  waves 
stopped  breaking  fthe  old  equilibrium  with  clear  water  above  a  few  t 

centimetres  of  dense  suspension  was  reestablished  in  about  ten  wave 
periods  time. 

The  two  above  mentioned  situations  of  the  heavily  plunging  I 

wave  and  the  gently  spilling  one  are  the  extreemes  with  respect  to 
breaker  introduced  mixings  and  their  effects  on  the  sediment  distri¬ 
bution  in  general  and  on  the  boundary  layer  in  particular  axe  pro-  I 

bably  very  different. 

The  laboratory  measurements  shown  in  Figure 5.20  are  from  under 
non  breaking  waves  and  spilling  breakers  of  the  same  height.  The  I 

breaker  turbulence  has  obviously  changed  the  upper  part  of  the 
profile  drasticaly;  but  the  lower  part,  that  is  CQ  and  axe 
unchanged,  so  the  concentration  magnitude  and  the  near  bed  vertical  ^ 

length  scale  l are  the  same  under  spilling  breakers  as  under  on 
breaking  waves.  This  tendency  is  not  contradicted  by  the  field 
measurements  of  the  present  study.  See  Figure  5.10.  ^ 

Hone  of  the  field  experiments  referred  here  we  carried  out 
under  really  heavily  plunging  waves.  All  the  experiments  marked 
"P"  for  plunging  in  Figure  5.10  were  performed  on  flat  straight  ^ 

beaches  where  the  upwelling  of  sand  with  entrained  air  never  be¬ 
comes  as  well  developed  on  shallow  bars  for  example,  nevertheless 
there  is  a  tendency  for  the  "P" -experiments  to  fall  above  the 


average  in  Figure  5.10  indicating  soma  injection  of  breaker  tur¬ 
bulence  into  the  Boundary  layer. 

A  theoretical  treatment  of  sediment  dispersion  due  to  wave 
breaking  does  not  seem  possible  at  the  moment  although  some  in¬ 
formation  about  the  turbulence  structure  can  be  gained  from  such 
studies  as  Miller  (1976) ,  Pe re grime  and  Svendsen  (1979) ,  Madsen 
(1981)  and  Thompson  (1982) .  Neither  is  it  possible  to  derive 
precise  quantitative  information  concerning  the  upper  part  of  the 
concentration  profiles  on  the  basis  of  presently  available  data. 

A  detailed  quantitative  discription  would  require  a  comprehensive 
laboratory  study  where  such  things  as  break  point  and  breaker 
type  can  be  clearly  defined.  These  are  never  constant  nor  well 
defined  in  the  field.  It  will  also  be  necessary  to  start  with  very 
uniform  sediment,  so  that  the  affects  of  settling  velocity  varia¬ 
bility  are  eliminated  in  the  first  place. 
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Figure  5.20 1  Concentrations  measured  under  spilling  breakers  and  non 
breaking  waves  of  the  same  height,  from  Nielsen  (1979),  h  •  0.4m, 

T  •  1.7s,  H  -  0.19m,  d  -  0.082mm.  The  spilling  does  not  change  con¬ 
centration  magnitude  or  profile  slope  inside  the  boundary  layer 
(z  *  2cm) ;  but  it  changes  the  profile  completely  at  higher  elevations 


INFLUENCE  OF  CURRENTS  CN  CONCENTRATION  PROFILES 


In  none  of  the  field  experiments  of  the  present  study  did 
the  steady  current  influence  the  suspended  sediment  concentrations 
significantly.  That  is,  none  of  the  data  show  a  significant 
tendency  for  Cq  to  be  larger  or  the  distribution  to  be  different 
due  to  a  stronger  relative  steady  current. 

The  relative  strength  of  the  steady  current  is  defined  here 


VR  »  (S2  +  v2)/(aw) 2  (5.23) 

where  the  mean  current  components  u  and  v  are  shore  parallel  and 
shore  normal  respectively  and  u  is  generally  measured  0.2m  above 
the  bed  while  v  is  measured  at  0.3m. 

The  range  of  VR  in  the  experiments  carried  out  over  the 
rippled  bed  in  the  bar  trough  at  Eastern  Beach,  Cippsland  was  0.02 
to  0.19  with  the  steady  current  being  predominantly  shore  parallel. 

The  strongest  longshore  current  was  0.32m/s  in  Test  54,  measured  0.2Sm 
above  the  bed.  The  current  veolocity  at  the  surface  was  about  0.5m/ a 
which  is  a  very  considerable  current  to  work  in.  Nevertheless  the 
ripples  on  the  bed  were  apparently  not  affected  by  the  current  at  all. 
They  were  sharp  crested,  shore  parallel  and  symmetrical ,  just  like 
under  a  pure  wave  motion.  This  shows  that  it  takes  a  very  strong 
longshore  current  to  change  the  bed  form  geometry,  the  boundary  layer 
structure  or  the  sediment  entrainment  process  over  a  rippled  bed. 

The  wave  boundary  layer  and  bed  forms  seem  to  be  somewhat 


more  sensitive  to  shore  normal  currents  like  rip  currents.  From 
visual  observations  it  seems  that  even  a  weak  rip  current 
(V  s  0.05)  tends  to  change  the  bed  forms  into  asymmetrical  dunes 
with  somewhat  larger  length  than  vortex  ripples  and  with  no  vortex 
shedding  from  the  landward  slope.  Unforunately  we  did  not  suceed 
in  getting  any  concentration  measurements  from  rip  currents.  This 
is  of  coarse  due  to  the  fact  that  rip  currents  are  very  difficult 
to  work  in.  The  laboratory  studies  of  Srevik  and  Aas  (1980)  and 
Kemp  and  Simons  (1982)  have  revealed  important  details  about  the 
boundary  layer  under  waves  and  a  strong  "shore  normal"  current 
over  a  solid  rippled  bed.  But  the  response  of  a  loose  sand  bed 
to  such  flows  is  definitely  worth  a  study  as  well. 
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SHORENORMAL  SEDIMENT  TRANSPORT 

Prediction  of  sediment  transport  under  waves  and  a  current 
is  a  very  complicated  matter  and  at  present  there  is  not  model 
in  existence  which  is  both  theoretically  based  and  in  agreement 
with  experimental  evidence. 

As  early  as  1963,  Inman  and  Bowen  showed  that  the  net 
sediment  transport  under  waves  and  a  codirectional  current  will 
just  as  often  be  against  the  net  current  as  with.  it. 

This  means  that  no  transport  model  by  which  the  net  trans¬ 
port  has  the  direction  of  the  net  flow  will  be  in  general  agreement 
with  the  empirical  evidence. 

More  such  evidence  has  since  been  obtained  by  many  detailed 
flume  experiments  at  the  Technical  university  in  Delft;  see  van  de 
Graff  (19801  All  of  these  experiments  as  well  as  those  of  Bowen 
and  Inman  were  carried  out  over  rippled  beds,  with  relatively 
weak  currents. 

The  explanation  for  the  tendency  of  the  net  transport  direction 
to  be  opposite  to  that  of  the  strongest  instantaneous  current  lies 
in  the  way  sand  is  entrained  over  ripples  (see  the  illustrations 
of  Bijker  at.  al.  1976) :  The  sand  that  is  picked  up  by  the  velocity 
in  say,  the  shoreward  direction  is  trapped  in  the  lee  vortex  and 
thus  not  effectively  moved  until  the  vortex  is  released  by  the 
flow  reversal  and  carried  seaward  by  the  following  seaward  current. 

In  this  way  the  weaker  offshore  current  carries  larger  con¬ 
centrations  which  is  likely  to  result  in  net  seaward  transport 
because  concentrations  grow  rapidly  with  the  entraining  velocity 
(like  u1*  or  u*).  The  picture  is  of  course  different  when  the 


current  becomes  so  strong  that  the  flow  never  reverses  through 
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most  of  the  water  column.  Further  complication  is  added  by 
the  fact  that  waves  and  steady  currents  have  very  different 
abilities  to  penetrate  the  lower  layers  where  the  concentrations 
are  largest.  Recent  laboratory  studies  by  Brevik  and  Aas  (1980) 
and  Kemp  and  Simons  (1982)  give  good  details  about  the  current 
distributions  under  waves,  over  solid  rippled  beds,  however  it 
must  be  borne  in  mind  here  that  the  bed  forms  in  for  example  rip 
currents  are  very  different  from  ripples  and  this  may  change 
the  turbulence  structure  and  the  current  distribution. 

For  quantitative  description  of  the  sediment  transport  it 
is  appropriate  to  write  both  velocities  and  concentrations  as 
Fourier  series 


u(z,t)  «  He{?  0  (z)  e1"**} 
s  m 


(5.24) 


c(z,t) 


itf 


C  (z)  a 
a  m 


imut 


} 


(5.25) 


since  the  whole  process  must  be  periodic.  In  general  u  (z)  and 

2S 

C  (z)  are  complex  numbers  of  which  the  argument  determines  the 
phase  shift  relative  to  the  fundamental  mode  of  the  waves. 

The  net  sediment  flux  through  a  vertical  element  dz  at 
elevation  z  is  then 


SQ^CzJdz  ■  E|o  {*) |  |Cn(z) j  cos  $  dz 


(5.26) 


where  $  is  the  local  phase  shift  between  U  and  C  given  by 

m  m 
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(5.27) 


»  Arg  (u  (z) }  -  Arg{c  (z) } 

is  is 

13113  dependence  of  the  sediment  'flux  on  local  phase  relations  is 
a  central  problem  in  non  steady  sediment  transport  modelling  and 
as  yet  there  is  not  model  which  is  capable  of  predicting  these 
phase  relations  on  theoretical  grounds. 

The  model  presented  by  Nielsen  et  al.  (1978)  and  Nielsen  (1979) 
predicts  these  phase  relations  as  well  as  net  sediment  fluxes  that 
aure  in  agreement  with  the  Dutch  measurements ,  with  appropriate  choice 
of  parameters.  However,  this  model  is  based  on  the  diffusion 
equation,  and  as  we  found  earlier  in  this  section  the  entrainment 
process  especially  over  ripples  is  not  a  diffusion  process. 


SHORE  PARALLEL  SEDIMENT  TRANSPORT 


In  the  direction  perpendicular  to  the  wave  notion  the 
problem  is  not  so  complex  because  there  are  no  oscillatory 
velocity  components  {except  for  maybe  pulsation  of  longshore 
currents  at  surf  beat  frequency)  so  the  total  sediment  flux 
is  simply 

h 

Q  -  /  v(z)  c(z)  dz  (5.28) 

o 

The  evaluation  of  c(z)  and  v(z)  is  not  straight  forward 
however.  The  concentration ,  c(z)  will  in  general  be  a  function 
of  u  and  v  as  well  as  of  caa  although  the  dependence  on  u  and  v 
seems  to  be  rather  weak  and  is  in  fact  hidden  in  the  scatter  of 
the  presently  available  field  data. 

Our  ability  of  predict  the  longshore  current  distribution 
is  also  very  limited.  Two  rather  recent  models  for  wave  current 
boundary  layers  are  those  of  Grant  and  Madsen  (1979)  and 
Christoffersen  (1982).  The  first  model  assumes  that  the  eddy 
viscosity  grows  linearly  both  inside  and  outside  the  wave  boundary 
layer  while  the  latter  assumes  a  constant  ^  in  the  boundary 
layer.  The  second  model  is  much  simpler  to  work  with  and  in  fact 
the  measurements  of  Kemp  and  Simons  (1982)  show  that  under  their 
experimental  conditions  is  more  like  a  constant  than  anything 
else  inside  the  wavs  boundary  layer. 

None  of  the  two  models  considers  the  possibly  vast  im¬ 
portance  of  the  orientation  of  the  vortices  over  rippled  beds. 

This  orientation  may  however  lead  to  very  different  eddy  viscosities 


felt  by  flow  in  the  shore  parallel  and  the  shore  normal  direction. 
Strong  vortices  with  3hore  parallel  axes  will  probably  cause 
much  more  resistance  to  longshore  currents  than  to  shore  normal 
ones. 

'Zhe  general  behaviour  of  bed  forms  under  combinations 
of  waves  and  currents  is  also  virtually  unknown  so  far.  We 
only  know  that  they  look  very  different  under  rip  currents  and 
longshore  currents  of  the  same  strength  relative  to  au. 
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